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ABSTRACT 

This  paper  develops  two  different  forms  of  asymptotic  series  for 
the  spatial  "behavior  of  steady  state  electromagnetic  fields.  The  "basic  vari- 
a"ble  in  each  series  is  the  wave-length  so  that  the  first  fevr  terms  of  each 
series  may  serve  as  excellent  approximations  to  the  steady  state  for  small  h. 

Each  asyjnptotic  expression  for  the  steady  state  field  is  obtained 
in  terms  of  what  is  called  the  pulse  solution  of  Maxwell's  equations.  This  is 
the  field  which  is  established  by  the  sudden  introduction  of  charge  into  the 
medium  and  which,  with  increasing  time,  approaches  therefore  the  static  field 
surrounding  that  charge.  Actually  the  discontinuities  of  this  field  and  of 
its  derivatives  with  respect  to  the  time  variable  determine  the  essential  part 
of  each  coefficient  of  the  tenas  in  the  asymptotic  series. 

The  result  permits  linlcing  the  asymptotic  series  with  the  field  fur- 
nished by  geometrical  optics.   It  is  shown  that  the  discontinuity  in  E  and  H 
which  exists  on  any  wave-front  created  by  the  sudden  introduction  of  a  charge 
or  by  reflection  and  refraction  is  precisely  the  value  of  B  and  H  given  by 
geometrical  optics.   It  is  shown  also  that  the  discontinuities  in  the  E  and  H 
o.f  the  pulse  solution  at  any  fixed  point  in  space  at  which  €  ^  ^  ^    and.  cr  are 
continuous  and  at  which  there  are  no  charges  are  discontinuities  on  wave  fronts. 
Hence  these  discontinuities,  which  are  given  by  geometrical  optics,  furnish  the 
amplitude  of  first  term  of  the  asymptotic  expansion  of  the  steady  state  field. 
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1.  Introdaction* 

This  paper  presents  two  forms  of  asymptotic  expressions  for  the 
steady  state  electromagnetic  field,  arising  from  a  periodic  source,  in  the 
form  of  series  whose  "basic  variable  is  the  wave-length.  The  series  need  not 
"be  convergent  Taat  can  at  least  furnish  an  asymptotic  representation  of  the 
field  which  is  useful  especially  at  small  wave-lengths. 

The  full  asymptotic  representations  are  huilt  up  from  a  knowledge 
of  what  is  called  the  pulse  solution  of  Maxwell's  equations.   This  solution 
is  the  field  established  hy  a  source,  i.e.,  charges,  having  the  same  geonetri- 
cal  distriljution  as  the  source  of  the  steady  state  field  being  represented  but 
having,   instead  of  a  periodic  time  behavior,  the  behavior  of  being  0  for  t  <  0 
and  1  for  t  _  0.   In  the  first  of  the  asymptotic  representations  the  first  term 
is  given  essentially  by  the  sum  of  the  discontinuities  in  time  of  the  pulse  solu- 
tion;  the  second  term  is  given  essentially  by  the  sum  of  the  discontinuities  in 
time  of  the  time  derivatives  of  the  pulse  solution;   the  third  term  by  the  sum 
of  the  discontinuities  in  the  second  time  derivatives  of  the  pulse  solution;   etc. 
Of  course  each  term  of  the  representation  contains  the  appropriate  power  of  the 
wave-length,  ^.  This  result  is  embodied  in  formula  (9.18).  Actually  formula 
(9.18)  gives  an  exact  representation  of  the  steady  state  field,  but  it  leads 
directly  to  the  asymptotic  representation. 

The  second  asymptotic  representation,  which  replaces  formula  (9.18) 
when  the  discontinuities  in  the  time  derivatives  of  the  pulse  solution  become 
infinite  is  obtained  through  the  use  of  a  series  of  complex  integrals  [formula 
(10.12)  J  .  Each  integral  depends  upon  the  singularities  of  the  integrand  which 
in  turn  depends  upon  the  discontinuities  of  the  pulse  solution.   Evaluation  of 
these  integrals  leads  to  series  in  the  variable  X, 

Both  asymptotic  representations  call  for  a  knowledge  of  the  discon- 
tinuities in  the  pulse  solution  and  its  derivatives.  For  the  purposes  of  this 
paper  it  is  assumed  that  the  f\Tll  pulse  solution  is  known.   This  solution  can  be 
obtained,  at  least  theoretically,  by  solution  of  Maxwell's  equations.  In  an 
application  of  this  theory  it  will  be  shown  how  the  pulse  solution  is  obtained 
in  another  manner  for  a  class  of  diffraction  problems, 

*  This  paper  was  written  up  for  Dr.  R.  K.  Luneberg  by  Morris  Kline.  Crudities  in 
the  exposition  are  therefore  to  be  charged  to  the  writer. 

l.This  application  is  now  being  written  up  as  a  report  to  the  Geophysical  Research 
Directorate  of  the  Cambridge  Field  Station,  U.  S.  Air  Force. 
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The  physical  reason  for  the  discontinuities  in  the  pulse  solution 
and  its  derivatives,  at  a  point  (x,y,z)  in  the  case  where  ^ ,  u ,  and  O"  are 
continuous  there  and  no  sources  are  present,  is  the  sudden  passage  through 
the  point  of  the  signal  as  it  comes  directly  from  the  source  and  the  sutsequent 
passage  of  other  signals  which  are  refracted  by  obstacles  elsewhere  in  the  medium. 

The  subject  of  this  paper  must  be  distinguished  from  mimerous  special 
studies  in  electromagnetic  theory  wherein  a  given  expression,  representing  the 
solution  of  a  particular  problem,  is  evaluated  approximately  by  obtaining  an 
asymptotic  form  of  this  expression  valid  for  small  wave  lengths.  For  exeur^jle, 
the  diffraction  of  spherical  waves  by  a  half -plane  has  been  so  treated.   The 
asymptotic  expression  obtained  is  limited  to  the  particular  problem. 

The  general  result  of  this  paper  makes  possible  a  clarification  of 
the  relationship  of  geometrical  optics  to  electromagnetic  theory.-'  This  rela- 
tionship is  established  by  means  of  a  theoretical  approach  to  geometrical  optics 

k 
which  was  developed  by  the  author  prior  to  this  time.   The  essence  of  this 

approach, described  as  the  theory  of  wave-fronts,  will  be  presented  here  for  the 
sake  of  continuity  in  presentation  and  for  convenient  reference.   In  brief  it 
says  that  the  geometrical  optics  values  of  zhe   electromagnetic  vectors  E  and  H 
at  £iny  point  x,y,z  axe   the  values  on  the  wave-front  of  the  field  created  by  the 
source,  or,  equally  well,  by  the  wave-front  of  the  pulse  solution  arising  from 
the  source.   Since  the  surfaces  on  which  the  discontinuities  in  the  pulse  solu- 
tion are,  at  points  where  £  •  /*-f  ^^^   ^   are  continuous  and  no  charges  are 
present,  wave-fronts,  it  is  apparent  that  the  first  term  in  the  asymptotic  ex- 
pressions is  essentially  that  given  by  geometrical  optics. 

The  usefulness  of  the  asymptotic  theory  of  this  paper  for  problems 
involving  small  wave  lengths  is  apparent.  The  use  of  even  the  first  two  terms 
of  the  asymptotic  expansion  will  furnish  a  better  approximation  than  that  fur- 
nished by  geometrical  optics  alone.   It  is  hoped  that  the  theory  herein  pre- 
sented will  be  applicable  to  many  problems  for  which  the  steady  state  solution 
is  either  difficalt  of  impossible  to  obtain  exactly. 


2.  See,  for  example.  Baker,  B.B.  &  Copson,  B.T.:  The  Mathematical  Theory  of 
Huygens'  Principle,  Oxford  Univ.  Press,  1939,  pp.  95-9S. 

3.  A  limited  discussion  of  this  relationship  may  be  found  in  Baker  &  Copson,  loc. 
cit.,p.  79.  See  also  Bom,  M.:   Optik,  Julius  Springer,  Berlin,  I933,  PP.U5-U7. 

U.  Luneberg,  E.K, :  The  Mathematical  Theory  of  Optics,  Brown  University  Lectures, 
19^,  published  in  mimeographed  form  by  Brown  University,  Providence,  R.I.   The 
germ  of  this  theory  of  wave-fronts  is  also  in  Born, loc. cit. See  also  Silver, 
Samael:  Microwave  Antenna  Theory  and  Design,  McGraw-Hill  Book  Co.  N.Y  I9U9 
Chapter  \,  .    f  ^  >. 
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2.   Basic  Facts  about  the  .Electromagnetic  Field. 

2.1  For  isotropic  media  Maxwell's  equations  are  written 

curl  H--    E^-     —    <rB=0 
(2.1)  °       '         ° 

curl  E  +-^  H^     =0 
c       t 

where  £  , ^  ,  cr  are  given  functions  of  x,  y,  end  z.  We  assume  that  these 
electromagnetic  parameters  are  sectionally  smooth  functions  of  x,  y.and  z 
that  is  that  discontinuities  in  the  functions  and  their  first  derivatives 
are  finite  jumps. 

When  sources  are  present  in  the  field  it  is  desirable  to  write 
these  equations  in  slightly  more  general  form.  Let  us  assume  that  electric 
currents  are  produced  by  an  enforced  motion  of  electric  charge.  This  means 
mathematically  that  the  current  field  consists  of  two  parts,  the  induced 
part,  which  is  given  by  crE,  and  the  enforced  part.  The  latter  is  the  cur- 
rent due  directly  to  the  motion  of  the  charges  of  the  source,  which  may  be 
at  a  point  or  distributed  in  space.      Vfe  represent  the  enforced  current  by 

j=- F   (x,y,z,t)  where  F(x,y,z,t)  is  a  given  vector  field  determined  by  the 

source  distribution.  It  will  generally  be  assumed  that  the  components  of  F 
are  sectionally  smooth.  When  the  source  is  taken  into  account  we  must  con- 
aider  instead  of  the  homogeneous  equations   (2.1)   the   inhomogeneous  equations 


(2.2) 


curl  H     -     -E^-~o-E     =     i-F, 
etc  c     t 

curl  E     +    -^  H^     =  0  . 
c     t 


The  problem  we  shall  consider,  which  we  shall  refer  to  as  the 
general  radiation  problem,  is  to  find  solutions  E(x,y,z,t)  and  E(x,y,z,t) 
which,  at  t  =  0,  assume  given  initial  values  E(x,y,z,0);  H(x,y,z,0). 

In  most  applications  we  shall  assume  that  F,  E  and  H  are  iden- 
tically zero  for  t  i  0.  This  corresponds  to  the  physical  situation  that 
an  electromagnetic  source  begins  to  operate  at  t  =  0  and  an  electromagnetic 
field  is  established  under  its  influence.   If  moreover  the  enforced  currents 
are  confined  to  a  small  neighborhood  of  the  origin  (o,o,o),  then  F  =  0  except 
at  this  point  and  the  resulting  electromagnetic  fields  represent  the  radiation 
from  a  single  point  source  at  x  =  y  =  z  =  0. 
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In  case  (r=  0  it  follows  from  the  first  of  equations  (2.2),  "by 
taking  the  divergence  of  "both  sides,  and  integrating,  that 

(2.3)  div  e  E  =  -  div  F. 

Since  div  6  2  =  ^  ,  the  charge  distriljution,  we  see  that  f   also  determines 
the  charge  distribution. 

2.2.  For  many  practical  applications  we  can  assume  that  the  enforced 
current  is  given  "by  a  vector  field  F(x,y,z,t)  which  is  periodic  in  time.   We 
expect,  for  physical  reasons,  that  an  electromagnetic  field  is  estahlished 
with  increasing  time  which  is  likewise  periodic  in  time  and  has  the  same 
period  as  the  source.  This  limit  solution  is  called  a  steady  state  solution 
of  Maxwell's  eq-uations.   (Its  estahlishment  will  be  demonstrated  mathematically 
from  the  solution  of  the  ahove  general  radition  problem.)  Since  any  time 
periodic  solution  can  he  obtained  hy  superposition  of  simple  harmonic  solutions 
we  may  assume  without  loss  of  generality  that 

(2.U)       F  ='fc(x,y,z)  cos  cut  +  t*  (x,y,z)  sin  cot 

where     c     and  c*  are  vectors   independent  of  t. 

The  associated  steady  state  solutions   then  have  an  analogous   form 

E  =  a(x,y,z)  cos  cot  •••  a*(x,y,z)   sincot 


(2.5) 


H  =  ^(x.y.z)  cos  cot  +  to*(x,y,z)  sinCOt 


where  a,    a*,  b,  b*  are  vectors   independent  of  t. 

The  discussion  of  steady   state   solutions   is  greatly  facilitated 
by  interpreting  the  vectors  F,   E,  H  as  real  parts  of  complex  vectors.      In  fact 
we  may  write 

F=Rf|e-^-*} 

(2.6)  E  =  Efue-^^*} 


H  =  r{ 


-».     -icotl 
v  e  ] 


where  g,  u,  v  are  complex  vectors 


g  -  c  +  i  c* 
(2.7)  u  =  a  +  i  a* 
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u  and  V  are  called  the  amplitude  vectors  of  the  electromag:netic   field  E,H; 
T  =  ^     the  time  period,     ^  =  ^  =     ^     the  frequency.  X  =     ^     the  wave 
length. 

V/e  ohtain  conditions   for   the  complex  vectors  u,  v  "by  introducing 
into  Maxwell's  equations   the  complex  vectors  whose  real  parts  determine  F 
E,   and  H  loy  (2.6);      this  sutstitution  yields   the  time  free  form  of  Maxwell's 
equations 

curl  "v     +ik€*u=     -  ikg 
(2.8) 

curl  u  -  ik^  V     =0 

Iflli 
where  6  •  =  €  +     -^^^     is  called  the  complex  dielectric  constant  and  where 

k  =  fi^  =     ^ 
c  X       ' 

In  those  parts  of  the  space  which  are  not  occupied  "by  sources  we 
have  g  =  0.  At  such  points  we  find,  merely  "by  taking  the  divergence  of  each 
equation,   that  the  steady  state  solutions  n,  v  satisfy  the  conditions 

(2.9)  div  e*u  =div//v  =  0. 

2.3.  Using  the  usual  meanings  for  D  and  B,  we  introdace 
W  =  ^  E  •  D   ,  the  electric  energy  per  unit  volume 

W  =  g—  E  •  E   ,  the  magnetic  energy  per  unit  volume 

W  =  V  ■•"  W^     ,  the  electromagnetic  energy  per  unit  volume 

S  =  ^  E  X  H,  the  radiation  vector  of  Poynting. 

The  time  period  of  electromagnetic  fields  is  in  general  so  small 
that  it  is  impossible  to  observe  the  actual  oscillations  of  B  and  H.  The  same 
rapid  oscillation  is  true  for  the  energy  density  W  and  the  flux  vector  S.  We 
are  hov/ever  in  the  position  to  observe  the  average  values 

(2.10)  ^  '     m  J      (  ^  ^^  *  >"  S^)  <1*  a°^  S  =  y|f^  y   (BxH)dt, 

o  o 

where  T  is  the  period  of  oscillation. 
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Let  us  express  these  observable  average  values  in  terms  of  the 
complex  vectors    ^     and  '^     .      We  introduce 

E=^(ue  +ue         ) 

(2.11)  ^ 

H=^(ve  +ve  ) 

in  formulas   (2,10)  and  find  readily  that 

W     =     5^(6/u/^+/i/v/2   ) 


(2.12) 


S       =      r^   (U   X   V    "♦•  U   X  v) 


where  u  is  the  conjugate  complex  vector  of  u,   that,  u  =  a  -  i  "S*,  |ul=V  a     +  %*    , 
and  u  X  V  =  (li  +  i  1*)  X  (t  -  i  S*)  =  a  X  f  +  i  a*  X  "5  -  i  a  X  S*  ■»■  a*  X  "5*  . 

2.U.    The  electric  vector 

(2.13)  E  =  S  cos  aPt  +  a*  sin  o:»t  =     |  (u  e"^'*^  +  u  e*^*^) 

descri.'bes  a  closed  curve  in  space  if  t  varies  from  0  to  T.   The  gdometrical 
characteristics  of  this  cujrve  provide  us  with  another  otserva'ble  characteristic 
of  the  electromagnetic  field,  its  polarization. 

We  note  first  that  the  curve  must  "be  in  the  plane  determined  "by  the 
vectors  a  and  a*.   It  thus  mast  be  normal  to  the  vector  a  x  a*  and  hence  normal  to 
the  vector  -rr-  (uXIT).   This  result  allows  us  to  determine  the  position  of  the 
plane  of  the  polarization  curve  directly  from  the  complex  vector  u. 

We  show  next  that  the  polarization  cuiTre  is  sn  ellipse.  From  equation 

(2.13)  it  follows  that 

-2i<^t  ^  -2   2icu»tx 
e      +  u-   e     ) 

-2ic*^t  ^  -2  2icotv 
e      +  u^  e     ) 

-2ioJt  ^  -2  2icotv 

+  u     e  )       , 

2iCL>t  -2icot 

These  are  three  linear  equations   m  the  two  unknowns  e  and  e  .     A 

common  solution  is  possible  only  if  the  determinaii.t  of  the  coefficients  and 


A- 

1 
2 

°1^      =     1. 

S^ 

A- 

1 
2 

-^^z   '   k 

^4 

A- 

1 

2 

Vj   =   5 

(u^  e 
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constant  terms  is  0  ;      that  is ,   if 


C? 


(2.U) 


1       '  2 

-     2^^1      »^ 

.2 


-2 


„2  1        -  2-2 

2  2     2  2     '2     •        2 


111.2  1       -  2  -2 


=     0. 


2  "*3"3    •  "^3     •     "3 

i.e.,    if  B  =   (E- ,E_,B   )   lies  on  a  surface  of  second  order.      Since  E  describes 

a  plane  curve  we  conclude  that  the  polarization  curve  mast  "be  a  conic  section. 

it  must  "be  an  ellipse  "because  it  is  closed. 

The  state  of  polarization   is  characterized  by  the  length  amd  the 

position  of  the  main  axes  of  the  ellipse.     We  find  these  axes  "by  the  extreme 

values  of  the  expression 

/-,  •..^^  ^2  1/2     -2i£ot  ^  -2     2ia;t     ^  ^         -v 

(2.15)  E      =     ^u     e  +u     e  +  2u  •  u)     , 

By  differentiation 


2     -2ia?t       -2     2icut 


and  hence 
(2.16) 


2iCc>t 


-     -V~2       • 
u 


The  extreme  vectors  E  then  may  "be  found  "by  inserting  the  value   (2.l6)  in  ex- 
pression (2,13).     The  lengths  of  these  vectors  follow  from  formula  (2,15), 
namely, 


(2.17) 
The  ratio 

(2.12) 


E^     =     |(u  .  u     +Vu%^     )     . 


^2     _     u  •  u     -  Vu    u" 


u  •  u    +Vu^  u^ 


determines  the  ellipticity  of  the  polarization.      The  polarization  is  circular 

2   -2 

1  =  u  =  0  ,   The  pc 

u  •  u  -  Vu^  u^  =  0 


if  K  =  1;  i.e.,  if  u^  =  u^  =  0  .  The  polarization  is  linear  if  K  =  0;  i.e.,  if 
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or  if  uXu     =0.      This   is  the  case  if  u  is  proportional  to  a  real  vector.     The 
ellipse  of  polarization,    then,    degenerates  into  a  section  of  a  straight  line. 

Since  Ti(z,y, z)  varies  with  x,y,z  we  see  that  the   state   of  polariza- 
tion will  in  general  "be   different  at   different  points   of  the  space.     A  wave  may 
"be  linearly  polarized  at  certain  points   in  space  and  elliptically  polarized  at 
other  points, 

3.   Discontinaity  Conditions  for   signal  and  Medium  Discontimiities. 

3.1.  Ordinarily  the  conditions  which  the  vectors  E,   D,   E,    and  B  must 
satisfy  at  a  surface  where   €  ,  ^  ,   or   cT   changes  ahruptly  are  derived  "by  a 
physical  argument  which   is  reasonably  satisfactory.     However,    these  conditions 
do  not  cover  the  case  of  discontinuities  existing  in  the   signal  or  field  itself 
at  points  where    ^,    ll,    ajid   cr  are  continuous.      Because  we  shall  he  interested 
in  the  mathematical  conditions  holding  on  a  wave-front,    that   is,    the  outermost 
region  of  penetration  of  a  signal,  we   shall  consider  signal  discontinuities  and 
in  doing  so  will  approach  the  whole  question  of  discontinuities  from  a  more  mathe- 
matical standpoint.      It  will  he  seen  that  these  considerations   lead  to  the  usual 
discontinaity  or  "boundary  conditions   arising  from  a  change  in  the   character  of  the 
medium. 

3.2.  VTe  derive  these  conditions  as  follov/s.      Consider  a  function 
A(x,y,z,t).     Let  G-  he  a  domain  in  the  four-dimensional  x,y,z,t  space,   and  let   the 

hyper-surface    1    he  its 
"boundary.     We  assume  that 
A(x,y,z,t)   is  continuous 
and  has   sectionally  contina- 
ous  derivatives.      Let 

^  =   ^^N'  ^N*    h'   S^  ^^  ^ 
four-dimensional  unit  vector"* 
which  represents  the  outside 
normal  on  1   ,     Then  we  have 
the  relation 


t 
i. 


■*■  % 


Fig.    1 


(3.1) 


/£..  E 


^ 


ds 


5,      If  the  equation  def  ining  F  is    ^(x,y,z,t)  =  0  then  x^=  ^0    ,y  =  X<P  , 
tjj  =  >N  ^ .    where  X  = 


"N 


^<P^' 


-V0x*<^y*   0z*0 
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where   dv  =  dx  dy  dz  dt  is  the  four- dimensional  volume  element  and  where  ds  is 
a  three-dimensional  element  of  the  hyper-surface  i     ,     The  analogous  relations 
wherein  x  may  "be  replaced  hy  y,z,   and  t  hold  also. 

The  proof  of  this  relationship  in  four-dimensional  geometry  is  merely 
a  generalization  of  the  familiar  fact  that   if  we  have  an  integral  of  any  f(x,y,z) 
over  a  volume,  we  may  integrate  partially  with  respect  to  z,    say,    substitute  the 
z-boundaries  in  the  integral  and  obtain  two  integrals  over  the  projection  of  the 
original  domain  onto  the  xy-plane.  However,   the  latter  two  integrals  gire  equiva- 
lent to  a  surface   integral  over  the  boundary  of  the  original  volume  in  which 
integral  the  integrand  must  be  multiplied  by  the   direction  cosine   of  the  normal 
to  the  surfgice. 

Consider  next  a  vector  field    X  (x,y,z,t).     We  find  that 

fjj]  curl  A  dv  =       Ilf  {U  xl  )  ds 

G  r 

(3.2)  ff/f  div  I     dv  =     J/f  (M  •  A  )  ds 

G  r 

////  X     6.r         =      ///   A     t         is 

G  *  r  ^ 

where  M  =   (x^.  J'jj.    Zjj.   0)  is  tl^e  projection  of  H   onto  the  x,y,z  space.      In  rela- 
tions   (3.2)  curl,    divergence,   and  cross  product  are  used  in  the  ordinary  three- 
dimensional  sense.     To  establish  the  relations  we  have  to  remember  merely  that 
curl  and  divergence  are  a  sum  of  first  order  derivative  terns  and  so  we  can 
apply  relation  (3.I)  to  each  term  and  then  recollect  terms  in  accordance  with 
the  meaning  of  cross-  and  dot  product. 

Let  us  now  assume  that  E,   H,   D,  B  and     B  are  continuous  and  have 
sectionally  continuous   derivatives.      If  we  denote  the  toteil  current,  namely 

0"  E  +  J—    F     by  I  we  may  write  Maxwell's  equations  as 

curlH.lD^=     ^(JE*1f,  ;   div  (I  *  ^  D^  )  =  0 


(3.3) 

ciirT   "Pi!  _        ^ 
C   ''t 


curl  E B.    =  0  :    div  B  =  0 
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We  integrate  these  equations  over  an  arbitrary  four-dimensional  domain  G 
and  make  use  of  the  relations  (3.2).   It  follows  that 


jj  j  (MxH  -  -I  D  -  i-  t„F)  ds 


c 


N 


=  ^  /#/ 


<rE  dv 


M 


(l*^D^)d8 


r 


B  ds  =  0  . 


r 

(3.^) 

yJJiHx^   +  -^  B)  ds  =  0 

r  ° 

These  integral  equations  represent  conditions  for  the  vectors  E,  H,  D.  B,   E 
and  F.  A  derivative  appears  only  in  the  second  equation  of  (3.^).   For  continu- 
ous fields  they  are  equivalent  to  the  differential  equations  (3.3)  insofar  as 
vfe  can  retrace  our  steps  and  derive  the  equations  (3.3)«   If"i  however,  the 
field  is  discontinuous  and  satisfies  the  integral  conditions  (3.U)  then  the 
differential  equations  (3.3)  will  still  he  satisfied  wherever  the  field  is 
continuous.  But  in  addition  to  this,  certain  conditions  can  he  obtained  "by 
(3»^)  which  refer  to  the  discontinuities  of  the  field  and  which  cannot  be 
derived  from  the  differential  equations.  Our  next  step  will  be  to  formulate 
these  discontinaity  conditions, 

3.3.  Let  us  assume  tliat  an  electromagnetic  field  v/hich  satisfies  the 
conditions  {j).k),   or  the  source,  is  discontinuous  on  a  hypersurface 


0(x,y,z,t)  =  0. 


We  choose  an  arbitrary  domain  G 
which  is  divided  into  two  parts 
G^  and  G  by  0=   0.  Also  the 

boundary  JT  of  G  is  divided  in 


two  parts 


^1  ^^ 


r_.  Let  r 


Fig.  2 


2*      -0 

be  the  part   of  0  =  0  which  lies  in- 
side G.      We  denote   the   discontinuous 
values  of  E,  H,   D,   B.CTE  and  F  on   P  "by 

»  ft 


E^,   H^,    D^,   2^,     crE^,   F^ 


and 


^2'  ^2' 


D^,   B^.     C^-E^,   F^ 


-12- 

aixd  assume  that  the  surface  normal  on    T 

o 

-— iriiiz:zz=iz==i=—     ^  0^*  0„.  <b„>   0+) 

V    0x  *  ^y  *  <^z  *  0t 

points   into  the  domain  Gp, 

We  consider  now  the  first  equation  (3.^)  aad  apply  it  to  the  total 
domain  G.      It  follows   that 

(3.5)         pV/(M^H-    ^    D-    I    t^)  ds  =    ^       ///fcris,    dv. 

However,    the  same  condition  mist  "be   satisfied  ty  *;he  domains  G,    and  G_  individu- 
ally.     This  gives,    using  grad  0  f or  the  three-dimensional  vector  {    0       0,0), 

dB 


///(MxH-^D     -|t^)     ds     .    .^(gradV^x  H^-  -%^.  I  Vr^^^,       ,       ,       , 
(3.6)  =    ^      ////, 


V0    +0    +0   +0. 
^x    ^y    ^2    '^t 


o-E     dv 


°i 


Tp  c  c  Tq  ^2cct2      /— 5 s s p 

■=    ^      IflUn    dv. 

We  add  the  equations   (3.6)  and  subtract  from  (3.5).      It  follows 

(3.7)  ///(grad  <p  x(H  -E    )  -  "^D  -D  )-  |   0,(r  -F   )J- ^°  =  0. 

r  c       1     .£       c       t     1     ^^/^2.  ^2.  ^2.  ^2 

However,    since    X      is  an  arhitraiy  hyper-area  of  0  =  0  we  conclade  that  the  inte- 
grand in  (3.7)  fflU-st  "be   zero  and  this  gives 

(3.8)  grad0x/H/-     ^  [d]  -     |    <P^  flj  -^     0 
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where     H     =  H^^  -  H^  and  D     =  D^^  -  D^  are  the  jtunp  discontinuities  of  H  and 
D  on  0=  0. 

We  apply  the  same  consideration  to  the  other  equatione    (3,U)   aixL 
thus   find  a  system  of  difference  equations; 

grad  0x  IE]   -     -—  [d/   -     1    (p^M   =  0  grad  JZ)   '  [ [l]  ♦     }^  h^j    =  0 

^^•'^  0 

grad  0  X  [e]  +  ~  [b]  =  0  grad  0  •  [b]  =  0 

which  establish  conditions  for  discontinuities  of  an  electromagnetic  field. 
They  are  the  counterpart  of  Maxwell's  differential  equations  end  can  he  oh- 
tained  from  them  hy  replacing  formally  the  differential  operators 

9       .S-       ^L         3- 
ar^  '  9y  •     az    '     St 

hy  multiplicative   operators 

A'  ^r  <^z'  *^t  • 

3.U.  A  very  special  case  of  these  discontinuity  relations,   which  will 
he  utilized  later  arises  when  the   source     1  is  discontinuous   at  t  =  0,    so 
that  the  hypersurface  0=  t  =  0  is  a  surface  of  discontinuilyof  F,E  aniH.It  follows 
from  the   first  two   of  equations   (3.9),    since  grad  0  is  now  identically  zero, 
that 

(3.10)  [CE  +       f]  =  0  and      fyWH?     =  0 
and  hence  since  B  =  H  =  0  for  t   <  0  that 

E(x,y,z.  +0)     =  -     I    F(x,y,z,   +  O) 

(3.11)  and 

H(x,y,z,   +  0)     =  0 

In  case   cr=  0  we  conclude  from  Maxwell's   equations    (3.3)  that 
•^  (diveB  •♦■     F)  =  0  at  any  point  where    £B  and  1  are  continuous   and  have 
continuous   derivatives.     From  (3.II)   it  follows   that   div  (es  +,  F)  approaches 
0  if  t  approaches  +  0  and  hence  div    es  =/«  =  -  div  F  in  a  certain  neighhorhood 
t   >  0.     However,   it  can  readily  he  shown  that  this   relation  /"  =  -   div  F  holds   for 
all  t   >  0  as  a  consequence  of  the  discontinuity  conditions    (3.9). 
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3.5.   The  discontinuity  relations   (3,9)  are  general  and  hold  regard- 
less of  the  reason  for  the  discontinuity  in  H,    D.   E,   B,   or  F.      We  shall     for 
example,    apply  them  to   the  surface  of  a  wave-front.      Discontinuous   electro- 
magnetic fields  appear  also  on  the  houndaries  of  material  "bodies.    In  this 
case  the  material  characteristics   of   the  medixim  are   discontinuous   functions 
on  surfaces  l^T  (x,y,z)  =  0  which  remain  fixed  in  the  three-dimensional  space. 
In  the  four-dimensional  x,y,z,t  such  a  surface  can  he  interpreted  as  a  cylindri- 
cal  hypersurface 

0(x,y,z,t)  =  ^(x.y.z)  =  0 
whose  generating  lines  are  parallel  to  the  t  axis.      From  (3,9)   it   follows  that 

grad>\r  X  [hJ     =0  grader  •[[!]+     ^  [dJ\  =  0 

grad\|^x  [i;]     =  0  gradxV  •     [  b]    =0 

Since  gradl^Tx  E  and  grad\^x  H  are   determined  "by  the  tangential  components   of  I 

and  H  and  gradlj/"*    (I  +     5^^^)  BJid  3  •  grad^^V  the  normal  components  of 

I  ■*■    ^  D^  and  B  we  can  express   (3.I2)  in  the  following  usual  manner, 

Thb  tangential  components  of  £  and  H  and  the  normal  compoitents  of 
I  ■••    5f  ^t  Q^'^  0^  ^  ^^  continuous  on  any  houndary  of  material  hodies. 

h.   Signal  Discontinuities. 

U.l.   This  next  article  will  consider  special  relations  vAiich  hold  on  a 
wave-front,    that  is,   on  the  outermost   surface  reached  by  a  signal  at  a  given 
time  t.      The  mathematical  laws   satisfied  on  a  wave-front  "by  the  surface  and  hy 
the  electromagnetic  vectors  on  the  front  constitute  geometrical   optics,   "ifhile 
this  concept  of  geometrical  optics   is  not  the  familiar  one,    it  will  he  seen 
that   the  laws  of  "behavior  of  fronts  and  associated  rays  (to  "be  considered  in 
the  next  article)  are  precisely  the  famllisir  laws   of  geometrical  optics.      Our 
approach  to  geometrical  optics  will  prove  useful   in  the  theory  of  steady  state 
solutions  to  "be  discussed,  later, 

k.2.   Consider  an  electromagnetic  field,    at   the  time  t  =  0,   which  is 

different  from  zero  only  in  a  certain  domain  G     bounded  "by  a  surface  T    .     The 

o  0 
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rig.  3 


electric  excitation  thus  shall  he  con- 
fined to  G  .   At  a  time  t  >  0  we  expect 
the  excitation  to  he  expanded  into  a 
large  domain  G,  hounded  hy  a  new  surface 
X  ,  Thus  to  any  time  t  there  helongs 
a  surface  4*(x,y,z,t)  =  0  which  determines 
the  houndary  of  the  electromagnetic  distur- 
bance. Such  surfaces  are  called  Wavefronts. 


A  set  of  wavefronts  given  hy  0(x,y,z,t)  =  0  for  different  values  of  t  in  the 
^,y,^-   space  generates  in  the  four-dimensional  x,y,z,t  space  a  hyper-surface.   Its 
"contour  manifolds"  t  =  const,   give  the  wavefronts  if  projected  onto  the  x,y,z 
space.   If  the  electromagnetic  vectors  E  and  H  are  different  from  zero  on  the 
originsd  surface  then  they  will  he  found  different  ^rom  zero  on  all  the  later 
wavefronts.  This  means  that  they  will  he  different  from  zero  on  the  hypersurface 

<P  =  0.      Since  the  field  is  zero  outside  this  surface  we  have  a  field  which  is 
discontimious  on  the  hypersurface  0=0.  The  discontinuities  [e] ,  f  HJ ,  f d/ ,  |bJ , 
[l]  moreover  are  equal  to  the  vectors  E,  H,  D,  B,  I  themeelves. 

In  this  article  we  shall  he  concerned  with  these  electromagnetic  vec- 
tors on  the  wavefronts,  i.e.  with  the  electromagnetic  field  on  the  hypersurface 

<P=  0.  Me   shall  find  that  these  discontinuities  estahlish  a  ^oup  of  phenomena 
v/hich  can  he  treated  separately  insofar  as  the  lav/s  of  propagation  are  independ- 
ent of  the  field  inside  the  wave  surfaces,  i.e.  inside  the  hypersurface  C^=  0, 
When  (^^  f  0  ve   may  solve  <|>(x,y,z,t)  =  0  for  t  in  the  form 

<^=4"U,y,z)  -  ct=0 
80  that  the  wavefront  at  the  time  t  is  given  hy  the  surface 

4^(x,y,z)  =  ct  . 
The  function l|/"(x,y,z)  is  called  the  wave  function.  The  electromagnetic  vectors 
on  the  wave  front  then  are  given  hy  the  vector  field, 

E*  (x.y.z)  =  E(x,y,z,  i  l^f ) 


H»  (x.y.z)  =H(x,y,z,  l\^)     . 
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E*  and  H*  thus  give  the   first   disturhance,    i.e.,    the  signal  which  is  observed 

at  a  point  x,y,z  at  the   time  t  =  — "Vf" (^»y»z). 

c 

Our  immediate  prohlem  will  he  to  determine  the  laws  satisfied  hy  E* 
and  H*  and  later  to  determine  them  from  some  initial  information. 

U.3.Let"^(x,y,z)  =  ct  he  a  set  of  wavefronts  obtained  hy  varying  t.     We 

characterize  the  velocity  v  of  propa- 
gation hy  the  velocity  of  a  point  P 
of  the  wavefront   on  the  associated 

orthogonal   trajectory,     SinceVK,\/r, 

X       y 

and  l//"    are  direction  numbers  of 
^  z 

V''=ctx  the  normal  to  the  surface  at  x,y,z,    the 

pj^       1^  velocity  vector  belonging  to  such  a 

point   is  given  by  the  vector 


(^•^)  (f'f'f  )     = 


"^l^vfWV, 


^M^x'^.V^z^ 


where  v  is  the  magnitude  of  the  velocity  of  propagation  and  wherein  x(t),y(t),  and 
and  z(t)  are  the  parametric  equation  of  the  orthogonal  trajectory.  Since,  on  the 
other  hand,  (see  Fig.  k) , 

YCP^)  -  -^iT^)   =  cCtg  -  t^) 

if  we  regard"^/(x,y,z)  as  a  function  of  the  x(t),  y(t)  and  z(t)  of  an  orthogonal 
trajectory  then 

dt~    ~    Yx    It     ■   ^y  dt"     '    Tz     dt 

dx 
dt 


dt  Tx     dt         ^y  dt  '  z 

By  substituting  for  — -  ,    etc.,   from  equation  (4,1)  we  obtain  the  relation 


(^.3)  VxKx*^l^y*y 


2  ^...2  ^     .2  c 


or 

(U.U) 


The  ratio 
(1^.5) 


V    = 


^/,    2    .  2~        2" 

V  TK^      Ty  *  Tz 


c 
n  =    — 

▼ 
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is  called  the  index  of  refraction  and  we  have  the  general  relation 


{k.S)  ^^"Wy"    V^z     '^^ 


2  a.       2  ^     .2     _     2 

Later  we  shall  find  that  n  is  determined  "hy  the  medium  alone,   independently  of 
the  particular  electromagnetic  excitation. 

The  vector 
(U.7)  ?    =     gradl^ 

is  called  the  wave  normal.   Its  components  are  denoted  hy 

p  =  -V^^.   q  =  A^y.   r  =  -g/^  . 

The  index  of  refraction  n  thus  is   determined  "by  the  absolute  value    Ip  I  of  the 
wave  normal     ^  , 

If  we  differentiate  the  wave  fiinction  \^/'(x,y,z)   in  a  direction  iO 
which  is  not  normal  to  the  wavefronts  we  ohtain  from 

ih.S)  ±W     -  -      ^^     +  .      dy  +  J.  —    -  c     —    -  ° 

d/o      ~  ■P    d,*'  cL^  i^  d/o  7 

the  wave  velocity  v  in  the  direction  p   ,  The  ratio 

f  =» 

is  called  the  index  of  refraction  in  the  direction  /<*  , 

It  is  given  "by  the  s cedar  product 
(U.9)  N=  p  .^ 

where  o  is  a  unit  vector  in  the  desired  direction  and  has  the  direction  cosines 

dx   d7   d  z 
d/<»  •  (^  •  d/*»   • 

U.U.The  electromagnetic  vectors  E*.  E*,  D*,  B*  on  the  wavefronts 

4>=    l/^Cx.y.z)  -  ct  =  0  give  directly  the  discontinuities  fB],  [eJ^    fDj,  [bJ 

of  the  field  on  the  hypersurface  ^=   0.  Consequently  the  equations  (3.9)  mast 

he  satisfied.  If  ?  is  continaous  on  such  a  front  then 


(U.10) 


0t 

grad  0  X  H*  -,  — -    D*  =  0 
c 

grad«^  X  E«  +  — -    C*  =  0  . 
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Since  4>  =  iK-  ct  it  follows  that 

grad  \^  X  E*  +  B*  =  0 
(H.ll) 

gradx//"  z  B*  -  B*  ^  0 

We  shall  first  assume  that  the  aeditua  itself  is  continuous,  i.e. 
that  the  quantities  ^  ,  /^  ,  and  <r  are  continuous  functions  of  x,y,z.  Later 
on  we  shaill  supplement  our  results  "by  allowing  discontinuities  edso  in  these 
functions. 

Under  this  assumption  ws  may  write 

(1+.12)  grad  -x/r  X  H*  *    6E*     =0 

grad  'iV  y.  B*  -  //E*     =  0 

We  conclude  readily  that 

(U,13)  E*  •  H*  =  E*  •  grad  y  =  H*  •  grad  "^  =  0. 

The  vectors  E*,  H*,  grad  lA"  from  an  orthogonal  system  of  vectors;  hence,  the 
vectors  S*  and  H*  are  tangent ietl  to  the  wavefront. 

The  equations  (U.12)  are  a 
homogeneous  system  of  linear  equations  for 
radT\^  the  components  of  B*  and  H"*.  Since  E*  and 

H*  are  assvuned  to  he  different  from  zero 
we  conclude  that  the  determinant  of  (U.12) 
must  he  zero.  This  implies  a  condition 
which  must  be  satisfied  "by  the  wavefront 
Fig,  5  function  "U/(z,y,z).  We  derive  this  condi- 

tion as  follows.  We  eliminate  H*  from  the  first  equation,  obtaining 

grad-y^  (gradl^x  E*)  *  e/^E*  =  0. 
grad^l/"  =  0  we  get 

(  e/^  -  grad^  V  )  E*  =  0. 


Since  E* 
(4.1U) 

and  hence 

(4.15) 


2       2     2 


y\fx  fy 

Wave  fronts  must  he  solutions  of  this  partial  differential  equation  which  is  known 
as  the  eiconal  equation  in  geometrical  optics. 

We  conclude,  with  the  aid  of  (H.  3): 

The  index  of  refraction,  n  =  —  ,  in  an  Isotropic  medium,  is  given 


by  the  relation 
(U.I6) 


n  =ye7i 


-iq- 

It  is  thus  independent  of  the  particular  propagation  and  also  independent  of 
the  direction  of  the  waves. 

U,5.  We  next  multiply  the  first  equation  (U.12)  with  E*  and  the  second 
with  H*.  We  find 

(U.I7)     (E*A  H»)  •  grad\|/  =  €  E*^  ,   (S*  X  E*)  •  grad'^/'  =  yU  H*^  , 

which  means  that  electric  and  magnetic  energies   are  equal  on  a  wave  front.     In 

formulae, 

6  E*^     =     l*nW 

(U.lg)  _ 

JUE*      =    Unw 

where  V  is  the  total  energy  density. 

Poynting's   radiation  vector  evaluated  on  the  wave  front,   namely, 

S*  =  ^  (E*A  E*)  must  be  proportional  to  gradl^   as  follows   from  (U.I3).     Prom 

(U.I7)  and  (U.lg)  it  follows  that 

S*     *   grad  y/  =  c\t 
and  thus 

(U.19)  S*  =     ^    V  gradl|/   . 

n 

This  gives  finally  "by  (4.15)  that 

(U.20)  |S*/  =  ^W  =  V  tf  , 

i.e.,    a  relation  between  flux  ajid  energy  density,   which  is  readily  Tinderstood 
physically. 

We  note  that  the  conductivity  CT  does  not  enter  in  our  results.     The 

shape  of  the  wave  fronts  thus  is   influenced,   also  in  metals,   only  by  e   and  x<  , 

5,   Wave-Fronts  and  Rays 

5,1. We  have  just  surveyed  the  principle  relationships  satisfied  by 
wave-fronts  themselves   as  well  as  by  the  electric  and  magnetic  vectors  on  tne 
wave-fronts.     We  shall  continue  the   study  of  wave-fronts  and  field  vectors  on 
the  fronts  by  introducing  a  family  of   rays.     A  ray  is  by  definition  a  curve 
which  has  at  each  point   in  space  the  direction  of  B*x  H*.      In  isotropic  media 
E*A  B*  has  the  direction  of  gradY''    where  \|/'(x,y,z)  =  ct  are  the  wave-fronts   (see 
formula  (4.13)).     Hence  in  isotropic  media  the  rays  are  orthogonal   trajectories 
of  the  wave-fronts. 
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After  obtaining  differential  equations  for  the  rays  we  shall  study 
the  "behavior  of  the  electric  and  magnetic  field  vectors  on  a  front  as  they  pro- 
ceed along  a  particular  ray  and  obtain  differential  equations  for  the  electric 
and  magnetic  vectors  along  rays. 

With  the  introduction  of  rays  it  "becomes  clearer  that  the  theory 
of  wave-fronts  is  equivalent  to  other  approaches  to  geometric  optics,  for  it 
is  possible  to  show  next  that  the  rays  satisfy  familiar  laws  of  optics  s\xch 
as  Fermat's  principle.   The  next  article  will  in  fact  establish  Snell's  lav/ 
of  reflection  and  refraction.  Moreover  it  is  possible  to  derive  the  mathema- 
tics of  wave-fronts  from  a  knowledge  of  the  rays.  Hence  wave-fronts  and  rays 
are  equivalent  approaches. 

We  may  repeat  that  the  essence  of  what  this  and  the  few  preceding 
articles  establish  is  that  the  discontinuity  relations  and  the  eiconal  equation 
axe  fundamental  and  that  given  a  family  of  surfaces  'S^  {x,y,z)   =  ct  satisfying 
the  eiconal  erjuation  and  two  vectors  E*  and  H*  satisfying  the  discontinuity  rela- 
tions (U.12)  on  this  family  we  have  geometrical  optics  studied  as  the  theory  of 
wave-fronts. 

5. 2. We  have  found  in  the  preceding  articles  that  electromagnetic  wave- 
fronts  "Vl/-  ct  =  0  must  satisfy  a  partial  differential  equation  of  first  order 

2     2     2  2 

which  for  isotropic  media  is  \i/  "*■  \l/v  *  V^   ~  ^^  ~  ^  '  This  equation  is 

called  the  characteristic  differentiad  equation  of  Maxwell's  equfi.tions.  The  con- 
structions of  wavefronts  thus  leads  to  the  ])roblem  of  integration  of  this  differ- 
ential equation.   It  is  known  from  mathematical  theory  that  the  integration  of 
partial  differential  equations  of  first  order  is  equivalent  to  the  integration 
of  a  certain  system  of  ordinary  differential  equations,  namely  the  characteristic 
differential  equations  of  the  first  order  partial  differential  equations.  Its 
solutions  are  certain  curves  in  space  which  are  called  the  bicharacteristics  of 
Maxwell's  equations,  or  electromagnetic  rays. 

We  introduce  these  curves  here,  not  by  mathematical  but  rather  by 
physical  considerations.  We  have  seen  in  article  k.k.    that  on  a  system 
'V|/(x,y,z)  =  ct  of  wavefronts  the  radiation  vector  S*  determines  a  vector  field 

•   —  ■ 

6.  See  the  Brown  University  Lectures  referred  to  in  footnote  3. 


(5.1)  i*  = 
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B*XH* 


in  space.  A  carve  which  at  any  of  its  points  has  the  direction  of  the  associated 
vector  s*  shall  be  called  an  electromagnetic  ray. 

Consider  a  set  of  wavefronts  ^  =  ct.     We  know  that  the  radiation 

vector  s*  is  normal  to  the  wavefronts. 
The  electromagnetic  rays  thas  are  the 
orthogonal  trajectories  of  the  set  of 
wavefronts,  i.e.  a  two  parameter  set 
of  space  carves  which  are  normal  to  the 
surfaces  \^=  ct.  We  represent  these 
p.   '  r      '  c-urves  in  parametric  form  x  =   x(  O") , 

y  =  y(  cr),  z  =  z( cr).   Since  they  are 
normal  to  the  surfaces  they  satisfy  the  differential  relations 

The  factor  X  can  "be  chosen  arbitrarily.  These  three  equations  form  a  system  of 
ordinary  differential  equations  for  x(cr),  y^O"),    z(  O").  Their  solutions  give  us 
certain  space  curves  which  are  geometrically  independent  of  the  choice  of  X.   In- 
deed, the  factor  X  influences  only  the  parametric  representation  of  the  curves  and 
not  the  geometrical  form.   If  another  system  of  wavefronts  ■^=  ct  is  considered, 
another  two  parameter  set  of  rays  is  obtained.   Thus  it  seems  that  any  solution 
^/ix,y,z)   of  the  eiconal  eq\iation  (U.I5)  characterises  a  set  of  rays  and  that 
the  totality  of  all  rays  in  a  given  medium  will  be  a  complicated  set  of  curves 
of  great  variety.  However,  we  shall  see  presently  that  this  is  not  the  case: 
The  total  set  of  rays  in  a  given  medium  is  determined  by  a  four  parameter  set  of 
space  curves  v;hich  can  be  characterized  independently  of  the  wave  functions 
^U(x.y.z). 

To  prove  this  significant  fact  we  divide  (5.2)  by  X  and  differentiate 
with  regard  to  CT ,   The  first  equation  yields 

do-  ^\d^^        ^xA  dcr     Tyx  ±<r  ^z/  dcr 
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The  other  equations  give  a  sljnilar  result  and  thus  a  system  of  second  order 

equations   is  formed: 

_d /  dx  s  ^      8  ,1     2v 

Xdo-  \icr'  9x  ^2  ^  '' 

d     I  d^z  \  _  _a_  /I   2v 


Ado-  Xd<r     3z  ^2 

in  which  the  particular  functionl//'(x,y,z)  no  longer  appears.  We  conclude  more- 
over from  (5.2)  that 

Hence  any  electromagnetic  ray  i(  cr),  y(  <r)  z( cr)  must  satisfy  the  differential 
equations  (5»3)  ^^^   *^®  condition  (5.^). 

Thewe  functions  x(  <r),  y(cr),  z(cr)  depend  on  the  choice  of  the  factor 
X.     However,  the  integral  curves  which  are  represented  "by  these  functions  do  not. 
Indeed,  a  curve  does  not  change  its  geometrical  shape  if  the  parameter  CT  is  su'b- 
mitted  to  an  arbitrary  transformation.   If  now,  on  a  particular  curve  for  a  cer- 
tain parameter  representation,  the  conditions  (5,3)  and  (5.U)  are  satisfied  then 
a  new  parameter  cr '  can  "be  introduced  such  that  the  conditions  (5.3)  and  (5.U) 
hold  with  X  =  1.  These  two  equations  are  the  "basic  eauations  for  the  rays  and 
could  "be  taiken  to  he  definitions  of  rays  in  an  introduction  to  geometric  optics 
which  is  independent  of  wave-fronts  and  "based  only  on  the  physical  properties  of 
the  medium  contained  in  the  quantity  n  =">^7^   . 

5, 3. Our  discussion  of  the  electromagnetic  vectors,  B*  and  H*  thus  far  has 
not  permitted  us  to  determine  them  but  Just  to  note  some  properties.  We  shall  now 
see  how  to  determine  the  vectors  E*  and  H*  on  a  set  of  progressing  wave-fronts 
'^=  ct. 

In  fact,  it  will  "be  seen  that  these  vectors  are  completely  determined 
along  a  particular  ray  if  they  are  known  at  one  particular  point  of  the  ray. 

The  vectors  E  =  E*  and  H  =  H*  on  the  wave-fronts  -^  =  ct  are  given  "by 
the  electromagnetic  vectors  E(x,y,z,t)  and  H(x,y,z,t)  if  we  insert  t  =  -  "^  (x,y,z). 
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From 

(5.5) 


it  follows  readily  that 


H*  =  H(x,y,z,  |a/^) 
7 

curl  £♦  =  curl  E  ■♦■  -    p  x  Ef 

c    *^         t 


(5.6)  1    ^ 

curl  H»  =  curl  H  +  -    p  x  H*     . 

c    •^         t 

Now  "by  Maxwell's  equations,  which  hold  everywhere  including  the  wave-fronts, 

curl  B  = UE. 

c  ~    t 

curl  H  =  -  (e  E.    +  1+n  crE). 
c  t 

We  replace  t  lay     —^  in  these  relations  and  then  suhstitute  in  (5.6).  Hence 

c 

one  ohtains  the  relations: 

€E»  +  p     X  H*     =  c(curl  H*  -     ^^L^  e*) 

(5.7)  *  * 

-J  X  E*  "♦■  /^H*     =  -c  curl  E*     . 

We  symmetrize  the  equations  (5.7)  "^V   adding  an  artificial  "magnetic 
conductivity"  k  so  that 
^5^gj  €E*  *P  xH*  =c(curlH»-  ^  E*) 

-p  X  E*  +  /«H*  =  -c(curl  E*  +  —  H^)  . 
X  V  c 

We  shall  of  course  make  k  =  0  afterwards.  But  the  symmetry  of  the  ahove  equa- 
tions allows  Us  to  derive  from  any  relation  a  second  one  simply  "by  replacing 
E*,  H*.  €  ,yU  ,  <r  ,  k  ty  H*.  E*.  -/i  ,  -€,   -k,  -cr  ,  respectively.  We  multiply 
the  first  equation  "by  /<  and  suhtract  the  vector  product  of  p  with  the  second  equa- 
tion.  It  follows  since  p  =  ^ /*  that 

-  (p  •  E*)  p  =  /«  curl  H»  +  I  X  curl  E*  +  —  [k(?  x  H*)  -  CT/U  B*7  . 
c      t  ,  0  ' 

or  since  "by  (4.7)  and  (U.12)  H*  =  ^  (p  x  E»)  and  p  x  H*  =  -  6  B»,  that  the  differ- 
ential expression  obtained  hy  suhsti tut ions  on  the  right  side,  namely, 

^5'9)     c^ri  (I  X  B*)  +  I  X  curl  E*  -  ~  (ke  +  o->*  )B* 
must  yield  a  vector  which  has  the  direction  of  j?. 

7.  In  the  equations  (5.6)  E^  means  partial  differential  of  E(x,y,s,t)  with  respect 


2k. 

c 
curl  E(x,y,z,t)  and  then  substitute  -^"^  for  t.  Likewise  for  curl  H. 


to  t  and  then  substitution  of  — t— -  for  t.  Curl  E  on.  the  right  side  means 

c 


Consider  now  the  vector  E*  x  curl     ^—     .      Since 

^       1  -♦         -►  1 

curl    p:  -  ^  curl  P     "*■  P     X  grad    ^     and  since  curl  f  =  0  by  virtue   of  f  =  gradV^ 

we  conclude  that 

E»     X  curl   ^    =  E*     x{3  X  grad  i)     . 

In  view  of  (U.I3),    the   right   side  equals   (E*  •   grad    rj)  p  and  so  the  left  side 
also  has   the  direction  of  ^.     Consequently  this   expression  may  he   added  to   (5,9) 
and  we  find  that 

(5.10)  curl  (E.    X  E»)  +     2.    X  curl  E*  ••■  E*  x  curl     2.  .     ifl     (k^      +   cru  )E» 

n  r*  M  CjU  ' 

has  the  direction  of  p. 

We  nov;  make  use  of  the  vector  identity 

(5.11)  curl   (AxB)  ♦  ^curl  E  +  Excurl  A 

=  -2(  A^  ^  *  A2  ^  +  A^  "fr^  "  ^  ^^"^  ^*^  ^^^  B+grad(A'E) 
which  holds  for  any  two  vectors  A  and  E.     We  introduce  A  =  2.     and  conclude,   since 
p   •   B*  =  0,    that 

(5.12)  2(p|f^    *q    IF    *^    |~)*E*>.div|    *    ^(ke    *    cr;u)E» 
must  have  the  direction  of  p  , 

Let  us  now  assume  that   the  rays  associated  with  a  set   of  wavefronte 
are  given  in  terms  of  a  parameter  Y  so  that  X  in  (5.^)  is  1.     Then  x{t  ) ,  y(  7*), 
2(7*)   satisfy  the  equation 

(5.13)  i^  +  J-^  +  z^  =  n^ 

and  hence  x  =  p,  ^  =  q,  z  =  r.  With  this  asaumption  we  may  write  (5.12)  as  follows, 
leaving  out  the  asterisk  for  themoment, 

(5.IU)         g;  *  I  E/*  div  I  +  |II  (k€   +  <r/^)  a  =  A? 

where  the  factor  X  is  as  yet  definite  hut  trndetermined. 

By  taking  the  scalar  product  of  p  with  the  equation  (5.1^)i  since 
p  •  E*  =  0,  we  ohtain 

^   2  ^p  dr     2  ^   dr 

n       *       n       * 
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But  since  p  =  i,  q  =  ^  and  r  =  z,  and  since  |p\  =   n  , 

^=p£ii.q|l+rfl=i  grad  n^  =  n  grad  n  . 
Hence 

(5.15)  A  =  -E  •  ^:4^  . 

We  finally  introduce  the  notation. 

The  operator  4y/V^T3ecomes  Laplace's   operator    ^'yf/' if  ^=  1. 

Substitution  in  (5.IU)  yields 

(5.16)         ||+    I  A^-/^E  +  U  •  ^^  )  gradx^   *    1^  (ke   *      <r/.)  B  =  0. 

A  similar  equation  will  te  found  for  H  "by  replacing  £  ,  li ,  CT  ,k  "by  -^  ,- ^  ,- k  ,-(77 
We  then  let  k  =  0  and  have  the  two  vector  equations  (restoring  the  asterisk) 

^  *  I  (A^T^   *  ^(T/W  )E*  *  (E*  .  i^  )  grad^  =  0 

^  *  I  (^e"^   *  ^Cr/^)H»  *  (H»  .  il^  )  grad-A-  =  0 


(5.17) 


wherein 

These  equations  are  differential  equations  along  a  given  ray  x(  f   ),y(  7*  ),z(  7*  ). 
They  allow  us  to  find  E*  and  H*  if  these  vectors  are  known  at  any  particular  point 
of  the  ray.  We  have  "but  to  suhstitute  the  quantities  6  ,  yU  ,  C"  t   n,  Zi  "*/^  •  ^°-<^ 
grad-vX^,  wherein  x,y,  and  z  have  the  values  x{Y),   y(  7* )  and  2(7*)  for  the  ray,  and 
then  solve  the  resulting  first  order  differential  equations  for  E*  and  H*. 

5,U.  We  can  simplify  the  equations  (5.I7)  considerahly.   It  follows  from 
the  definition  of  S,  and  equations  (U,19)  and  (5.I)  that,  in  an  isotropic  medium, 
p  =  €/i8*.      Consequently 

^i^"/^=  yU  div^  =  /^  div  €  t*  =  n^  div  ^*  +  s*. grade   . 
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Tbis  gives,  throughyUs*  ~  T    ^  and '^  •  gradC  =  -7-57  (since  the  rays  have  the 


direction  of  "f)   that 


d^ 


^i^f^=  n^  dlv  t*  +     ^   log  e 


and  similarly 


>^     dr 

^^=  n^  div  t*  +    Ip  log/^   . 
With  the  aid  of  these  expressions  we  write  (5«17)  as  follows: 

'^^"^J*^  +  i  (n^  div  t*  *  ^CU  )   (VT  E»)  *  (Ve  S*  .   £Iii^)gradV.  =  0 


(5.18) 


dr 


ii:^|!2  *  I  (n^  div  t*  *  ^O-JU)   (^  W)  *  (^  H»  •  i£^)  grad>^=  0. 
80  that  V6"S*  and  y^H*  satisfy  the  same  differential  equation. 

We  now  introduce  two  vectors  P  and  ^  by  the  equations 

VTe*    =  v^mw     P 

(5.19) 


^E*    =  VlffTW       Q 
where  W  is,    for  the  moment,  an  unknown  function.     These  vectors  have  the  same 
direction  as  E*  and  H*  Ifut  different  length.     We   expect  them  to  he  unit  vectors. 
We  introduce   (5.19)  in  (5. 18)  and  get 

1  P(i^    *  n2  divl-.  ^<r/.)  .    fL.*  (LJUE^l^)  g,,,^=  0 
and  an  identical  equation  for  Q.     We  now  submit  W  to  the  condition 

^5.20)  d         log  w  =  -(n2  div  t*  *  ^<r/4) 

and  thus  obtain  for  P  and  Q  the  simple  equations 


(5.21) 


dp 
dr 


^    P  •   grad  n  ^    1^       ,. 

+    s_ gradAf^  =  0 


We  shall  integrate  equation  (5«20)  after  establishing  a  helpful 
subsidiary  relationship.     Consider  a  "tube"  of  rays,  v;hich  as  indicated  in  Fig. 7, 

•  P  is  bounded  by  two  sections  w'   and  w  of  wavefronts. 

To  each  surface  element  dw  of  w  there  belongs  a 
corresponding  surface  element  dw'   of  w',   namely 
the   one  which  is  cut  out  by  the  rays  from  dw. 
Let 
*v  (5.22) 


K 


dw'  =     =T-    dw 


S-ig.   7 
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where  K(  T )   is  a  function,  vftiich  measures  the  expansion  of  the  wave   fronts  along 
e  ray.     Since  by  (5»i3)   'the  geometrical  line  element  d<r  on  a  ray  is  given  'by 
d  cr  =  ndT'     we  can  write   the  volume  element  dxdydz  in  the  form 
(5.23)  dxdydx  =  dw  dcr  =  ndwdr     . 

Now  let  us   apply  Gauss*    theorem  to  our  tube  of  rays   thus: 

JJJdir  t*  djcdydz  =    j]   t*  -  t«dw»   -      Jj    t*  '   n  dv  =    //  \T*\  dw'  -  // /s^/dw  , 

w'  w  w»  w 

where  n  and 'n'   are  vmit  normals  to  dw  and  dw'   respectively.      Since  ^=  n  e*  and 
\^1  =  n  ,  I  3*1  -     —     .      Hence,   with   (5.22),   the   last   expression  equals 

(         ^  -     _— i ^    kCT   )  dw     (where  T    is   the  veaue  at  w  and  T,   the 

value  at  w' )  'jj 


■llh 


d      t      1      \u  r   )r^     n*^^^*^ 


(^^^rr)  ^(  To)n 


Tr^skir)^        °  n2 

Hence  it  follows  from  (5.23)  that  at  "I*  (  T^) 
(5.2H) 


div  s*  = 

1-  ir^io6(-^))j 

f    i 

•^^=r^ 

With  this 

;   result 

(5. 24)  we  may  write 

(5.20) 

as 

nX  ~ 

W 
0 

UK      « 
0  0 

C     "^0     '^ 

dr 

(5.25) 


The  factor  W  fomally  introduced  in  (5.19)  is  the  electromagnetic 
energy,  for  from  (5.19 )  ^y   =   ,  ^»]Sfow  P  and  (J  can  he  made  to  "be  unit  vectors 
by  choice  of  the  constant  of  integration  in  (5. 20).   Since  '&*)(  H*  has  the 
direction  of  ^  ,  we  have  that   ij^   =  -3^-=_  =  -2-  =  t*    Reference  to  W 

in  article  2,3  shows  that  the  W  introdaced  here  is  electromagnetic  energy  if 
the  constant  of  integrntion  in  (5.20)  is  properly  chosen. 
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Several  other  implications  of  eopiations  (5.21)  \fill  be  noted.   We  form 
the  scalar  product  of  grad'^  with  the  first  of  equations  (5.21).   It  follows  that 

grad-^  •  ^  +  n(P  •  grad  n)  =  0. 

However,   n  grad  n  =     -r-^    grad"*^    ;   hence 

^5-^6^  1^     (P  .    gradV^  )  =  0. 

Thus,  if  (P  •  grad"^)  =  0  at  one  particular  point  of  the  ray,  then  it  will  "be 
zero  at  any  point  of  this  ray.   The  same  is  true  for  the  product  (Q  •  grad'^). 

We  next  multiply  the  first  of  equations  (5.21)  "by  P  and  use  that  P  •  grad-^s  0, 

dP  2 

It  follows  that  P  *  j;p-    =  0,   i.e.,  P  =  const,  along  the  ray.   Thus,  if  P  is  a 

unit  vector  at  one  point  of  the  ray  it  will  te  a  unit  vector  at  all  points. 

Finally  we  get  from  (5.21)  that 

^    dr       dr     dr      ^     • 
i.e.,  the  vectors  P  and  (^  are  normal  to  each  other  if  they  are  normal  at  one  par- 
ticular point  of  the  ray. 

From  (5.25)  it  follows  that  W  =  0  on  all  points  of  a  ray  if  W-  =  0 

W 

at  a  particular  point.     Also  if  the  medium  is  non-absorting  thenC=  0  pud    —^ 

is  constant  along  a  ray. 

Finally,  W(7'  )   is  completely  determined  "by  its  value  W(  7*  )   at  a  par- 
ticular point  along  a  ray.     We  can  thus  construct   the  vectors  E*  and  H*  on  any 
vfavefront")^  =  ct  if  they  are  given  on  one  particular  wave- fro nt '/^=  ct     .    . 

This  result  means  that    if  we  have  a  signal  such  that  discontinuities 
of  the  electromagnetic  field  are  different  from  0  only  on  a  sectionX^of  a  wave- 
front  then  the  rays  through  P  determine  the  region  where  the  signal     will  "be 
o"bserved.     This  does  not  exclude  the  possibility  that  the  signal  will  penetrate 
into  the  region  of  geometrical  shadow.     However  this   diffracted  excitation  does 
not  have  a  sudden  discontinuous  "beginning,    that   is,    it  is  not  given  "by  geometrical 
optics,   with  vdiich  we  are  now  dealing. 

6.   The  Behavior  of  Wave-Fronts  at  Refracting  Surfaces. 

6.1,   The     theory  of  wave-fronts  presented  thus  far  has  assumed  that  ^,/<, 
and  <r  are  continuous  and  have  continuous  derivatives  and  the  wave-fronts  have 
therefore  progressed  continuously  through  space.      This  article  v;ill  consider  the 
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effect  of  discontinuities  in  ^ ,  JU  t    and  Con   the  behavior  of  the  wave-fronts. 
Physically  the  presence  of  obstacles  or  of  lenses  in  a  medium  introduces  such 
discontinuities.   The  surfaces  on  which  the  discontinuities  in  ^  , // ,  and  (7* occur 
are  not  limited  to  planes.  Nevertheless  it  is  shown  that  the  usual  laws  of  re- 
flection and  refraction  are  obeyed  "by  the  wave-fronts.  Moreover,  the  reflection 
and  transmission  coefficients  are  obtained,  and  these  prove  to  be  identical  in 
form  with  the  Fresnel  formulas  for  reflection  and  transmission  of  a  plane  wave 
at  a  plane  surface. 

6.2.  The  effect  of  discontinuities  in  €  ,  ju  ^    and  CT   is  to  split  up  the 
wave-front  into  several  parts.   In  addition  to  a  transmitted  v/ave-front  there 
is  a  reflected  one.   If  several  refracting  surfaces  are  present  then  a  rather 
complicated  system  of  multiple  wave-fronts  which  overly  each  other  is  established 
in  time.     This  overlapping  is  avoided  if  the  wave-front s-)^=  ct  are  repre- 
sented in  the  four  dimensional  x,y,z,t  space  by  a  hyper-surface  <P --i^  ~   ct  =  0. 
It  consists  of  a  number  of  branches  as  is  schematically  illustrated  in  Fig.  8 
for  the  3c,y,t  space.  The  branches  are  Joined  together  in  groups  of  three  on 

the  refracting  surfaces,  however 
at  different  "heights",  i.e.,  at 
different  times. 

Consider  at  such  a  triple 
of  branches  joined  together  on  the 
refracting  surfaceJC*  : 

^  =  f(^   .7) 
(6.1)  y  =  g(^    ,.7  ) 

z  =  li(^    .9  ) 

The  three  hyper   surfaces  given  by: 
(6.2)U^  ct  =  0,y*-  ct  =  0,andT^'-  ct 
represent  the   incident,  reflected, 
and  refracted  wave-fronts,    res- 
pectively. 


=  0 


Fig.   9 
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In  order  that  these  hyper  surfaces  he  Joined  together  on  J*  we  nmst 
require  thatT^^*"/'*     ^  yf/-*  oar   .      This  Keans   that 

(6.3) 

If  we  differentiate  these  equations  with  respect   to  ^    or/7  we  ohtain  interior 
derivatives   of  ■^,  y^*,  l^'  which  again  raust  he  equal.      It  follows   that 


(6.U) 


and 


wave  normals 


^v-i  -y^x^^^   *  (^;  -y-y^«^  *  (v-i  -^z)^f  ==  0 

These  equations  show  that  the  differences  p*  -"p  and  p"'  -  "p  of  the 

are  parallel  to  the  normal  M  of  the  refracting  surface  at  the  point  «  ,  ^  .  We 
write 

p*=T+  r,  M 
P'  =  P  •»■   r^t 

and  conclude  that  all  four  vectors  M,  "p,  p"*,  p"'  lie  in  one  and  the  same  plane, 
the  plane  of  incidence.  We  also  conclude  from  (6.5)  that 

p*  X  i«  =  p  X  tr 


(6.6) 


p*     X  M  »=  p"  X  M 


By  taking  the  ah solute  value  of 
these  equatioiis   and  rememhering  that 
i2  _   ,-»..,  2  ?. 


2  ,2 

=  n'     we 


get  the  law  of  reflection  and  Snell's 
law  of  refraction: 

(6.7)       sin     Q*  =  Bind 
n'    sin     ^  •  =  n  sin  ^ 


-31- 


The  physically  consistent  solutions  of  these  equations  are 

B*  =     V  - ff 

and  an  angle  Q  '   of  the   interval  0     i    ^    =     tt/2     provided  that  Mis  oriented 
relative  to  p  such  that  0^0=       tt/2  ,     The  eoTiation  n'   sin   ^  '  =  n  sin^ 

has  no  solution  if  — r     sin^  >  1.      This   is  the  case  of  total   reflection  where 

n' 

the  signal  is  not  transmitted  through  the  refracting  surface. 

6,3»  We  study  next  the  electromagnetic  field  in  the  neighborhood  of  the 

manifold  J  where  the  three  hyper- 
surf  aces  i^=  ct  =  0,  •<//-*   -  ct  =  0 
and  -j^*  -  ct  =  0  intersect  each 
other.  This  manifold  J  lies  on  the 
cylindrical  hypersurface  Z^ {x,y,z) 
which  has  the  refracting  surface  as 
"basis.  The  neighborhood  of  J  is 

divided  into  five  regions  hy  the 
Fig.  11 


=  0 


*-  X 


four  hyper-surfaces.     Let  us  assume 

g 
that  the  vectors  E  and  H     are   con- 
tinuous  in  these  regions  and  that  they  attain  finite  limits 

^0'  \'  h'  ^3'  \ 

^o'  h'   ^2*  ^3*  ^1; 
if  point  P  of  the  manifold  J  is  approached. 


-^'=ct 


r=o 


Pig.  12 


We  denote  the  surface  nonaal 
of  the  refracting  surface  "by  M  and 
consider  the  limits  il  ,B  ,B  E^  and 
H^.H  ,H  ,H^^.  By  applying  the  con- 
ditions (3.12)  derived  in  article 
3.5  we  find  that 

(E^-E^)  X  K  =  0   (H^-H^)  X  M  =  0 
(6.8) 

(Bi^-E  )  X  M  =  0  (Hj^-H  )  X  M  =  0 


8.   We  shall   denote  the   incident  wave-front  "by  E  and  H  instead  of  E*  and  H*  and   shall 
use  E*  and  E*  for   the  reflected  wave-front   in  this  article. 
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We  next  consider  the  boundary  values  of  the  discontinuities  on  the  hype rsurf aces 
>/^=  ct,  -\^*  =  ct,   ->J^  =  ct.   They  are  given  "by 

0    4-  0    4 

(6.9)      E*=.\-\  H*  =  H,-H^ 

E«  =  Eg  -  E,  H'  =  H^  -  H 

and  represent  the  incident,  reflected  and  transmitted  signal  at  P  respectively. 
It  follows  readily  "by  (6.8)  that 

(E  +  E*)  X  M  =  E«  X  M 
(6.10) 

(H  +  H*)  X  M  =  H«  X  M 

These  discontinuities   satisfy  also  the  vector  equations 
^  X  H  +   6  E  =  0,      p»  X  H*  +  €  E*  =   0,  p'   X  H'   +    C  'E'  =  0 

(6.11) 

p  X  E  -/U  H  =  0,      p"*  X  H»  -/^H*  =   0,  p»   X  E'    -   /<'H'   =  0 

where  'JT  =  grad-^  ,  p*  =  grad  ■^*,  p'  =  grad  -^'  are  the  wave  normals  of  the  three 
wave-fronts. 

6.4,  We  have  seen  that  these  three  wave  normals  and  the  v  ?tor  M  lie  in  one 
and  the  same  plane  normal  to  a  unit  vector  S  defined  "by 

(6.12)  pxM  =  p*  xM='p''xM  =  kS 

The  constant  k  equals 

k  =  n  sin^      =nsin     0  *  ~  n'  sin  ^' 
In  addition  to  S  v^e  introduce  three  other  unit  vectors 

(6.13)  K  =  Sx^,K*     =Sx    ^     ,N'=Sx     ^ 
^    •    -"  n  n       •  n' 

which  are  normal  to  the  fncoming,   reflected  and  refracted  ray  respectively.     The 

vectors  U.S;   N*,   S;   K'.S     thus   form  three  pairs  of  orthogonal  unit  vectors  at 

right  angles  attached  to  each  of  the  three  rays.      Conseopiently  we  may  express 

the  vectors  E,H;      E*,H*;    E',H'   as  linear  combinations  of  the   associated  vectors 

N  and  S.     We  verify  readily  that  the  equations   (6.11)  are  satisfied  if  we  let 

V€B  =  «X  N  +  /9S    ;     VS"  E*  =   oc  ♦N*  +  ^  *S   ;     Vc'E«  =    <YiN«   +   /5  'S 
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V/e  shall  shortly  introduce  these   expressions  into( 6. 10).      We  have  first  from  (S.lU) 

that 

Ve  (E  X  M)  =o<(H  X  M)  +  ^(S  X  M) 
Since  by  relations  (6.12) 

N  X  M  =  (S  X  £•)  X  M  =  -  COB  ^  S 
we  get 

(E  X  M)=  ^^  cos  ^  S  +  ^j^  (S  X  M) 
and  similarly 

(6.15)  (H  X  M)  =  -^  cos  ^  S  +  -0^   (S  X  M) 

Corresponding  equations  are  found  with  regard  to  the  vectors  E*,H*  and  E',E'. 
Thus  equations  (6.10)  yield,  en  account  of  cos  0*  =  -   cos 5   ,  the  four  condi- 
tions   

which  can  "be  used  to  express  ^*,    ^ '»;  e<',  jff*   in  terms  of  the  incident  com- 
ponents 0(  and  ^  .   The  result  is 

^fiK  ^r^.  g°a  g'         ^fT'  ,/7^  cos  <$>' 

« *  _  V  A'  "  y^^  cos  S  S*     ^^~e^  ~  y~M*     cos  ^ 


(6.17) 


V/<»      y  £«    cos  ^  v^i    y^i    cog  ^ 


of 


^fK^JlLss±.Al.        ' '^       ./T ^JTT  cos^' 

y/<'  ¥€'003^  V   ^>       V^i      cos  ^ 

In  relt^tion  to  the  plane  of  incidence  it  is  customary  to  call  the 
components  of  the  electric  discontinuity  E  with  regard  to  N  and  S  the  parallel 
and  normal  component  of  B.  One  uses  the  notation 


(6.18)  R  =  -2L^  R  =  -^ 

D  =  -2 —  D  =  -S_ 


It   follows  letting//r/<,=l;  V£"=  n,     V?=  n'   that 


n       cos    Q '  n  cos    <?• 


(6.19) 


R 
p 

=  A 

P 

X 

ri'      cos   Q 

E«A„ 

B       S 

n" 

n 

n« 

■f 

cos  Q 

1 

■♦• 

n     cos   Q  • 
n*   cos  S 

cos    0  ' 
cos  0 

D 
P 

=  A 
P 

'h 

> 

D  =A, 
s     s 

n' 

1 

+ 

n       cos  S  ' 
n'      cos  ^ 

n 

+ 

cos    0  ' 
cos  0 

The  alDove 

formulae   are 

identical  with 

Fresnel's   fonmilae 

reflection  and  transmission  of  a  plane  wave  on  a  plane  surface  or  refraction. 
Their  significance,  however,  is  more  general  since  they  also  give  the  reflection 
and  transmission  of  any  discontinuous  signal  along  any  ray  which  strikes  a  sur-, 
face  of  refraction  of  arbitrary  shape.  We  have  to  consider  the  equations  (6. 17) 
as  the  supplement  of  the  transport  equations  for  signals  in  a  continuous  medium. 

We  remark  explicitly  that  these  formulae  are  valid  also  if  two  absorb- 
ing media  are  separated  by  a  surface  of  refraction.   The  conductivity,  C ,  thus  is 
found  to  be  without  effect  upon  the  reflection  and  transmission  on  the  surface  it- 
self.  It  plays  its  part,  however,  in  damping  the  signal  when  it  proceeds  from  the 
point  of  incidence  into  the  interior.  We  shaill  verify  this  result  later  by  exact 
integration  of  Maxv/ell's  equations.   Though  ^  has  a  considerable  influence  upon 
the  steady  state  solution  it  is  without  direct  effect  on  the  transient  phenomena 
on  a  proceeding  wave-front. 

7.  The  Asymptotic  Development  of  Steady  State  Dipole  Fields. 

7.1.  The  preceding  articles  dealt  with  the  first  excitation  which  reaches 
a  given  point  in  space  from  an  electromagnetic  source.  This  first  excitation  or 
wave-front  follows  certain  laws  which  are  quite  independent  of  the  field  estab- 
lished afterwards.   The  relation  of  this  theory  to  the  problem  of  the  paper  will 
soon  be  apparent. 

We  turn  now  to  the  steady  state  fields  set  up  by  the  source,  and  we 
shall  investigate  the  problem  set  for  this  paper,  namely,  to  obtain  asymptotic 
expressions  for  the  steady  state  fields  in  terras  of  the  variable  k  =  r —  , 

A 
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Before  developing  the  general  theory  we  shall  illustrate  the  nature 
of  it  "by  treating  the  field  of  a  dipole  in  a  homogeneous  as  well  as   in  a  non- 
homogeneous  medium.      The  mathaaatics   of  the  dipole  field  used  here  is  not  new. 
However,  we  shall  regard  it   in  the   light  of  our  general  theory  and  shall  gain 
insight   into  the  significance  of  the   terms. 

It  will  he   apparent  that  for  small  wave  lengths  the  amplitudes 
of  the   field  vectors   of   the  steady  state  field  satisfy  the  same   conditions  as 
do  the  electromagnetic  vectors   on  a  wave-front, 

7,2,     To  treat   steady-state  fields  we  return  to  Maxwell's   equations  as 
presented  in  article  2  end  to  the  various   formulations  of  these  equations  given  in 
that  article. 

We  shall  assume  first  that  we  are  dealing  with  a  homogeneous  medimn 
in  which,   for  simplicity,  6  -  M~  ^  a°^<i  ^-  0. 

Enforced  periodic  currents  shall  he  given  hy  the  vector  field 
(7.1)  ^    R(F^)  =  r(-   icof  e-^*-*^]      . 

where  F  =  g(x,y,z)e~^**'*.      The  steady  state  field  v/hich  is   estahlished  under  the 
influence  of  this  current  distrihution  must  satisfy  eq\iation  (2,6)  which   in  our 
special  case  reduces   to 

curl  V  +     iku  =  -ikg 

(7.2) 

curl  u  -  ikv  =  0, 
We  solve  these  equations  as  follows.   We  define  first  a  vector  fiela 
A(x,y,z)  "by  the  integral 

^7.3)         A^-^flfl-d  .^.^  )4^  d^  d^di; 
where  r  =V^(x-  ^.  )^  -^  (y  -tj)^   +  (z  _  i;  )2 

Let  us  assume  that  g(^     ,^  ,  4  )  is   sectionally  smooth  i.e.   that   any 
finite  domain  may  he  divided  into  finitely  many  suhdomains   in  which'g'is  continu- 
ous and  has  continuous   derivatives.      Then  in  any  of  these  suhdomains  ^satisfies 
the  differential  equation  in  x,y,   and  z, 

(7.U)  ^Tf  k^  T    =  ikl"  . 

This   statement  may  he     verified     hy  the  use  of  Greea's  theorem. 
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We  nov;  construct  two  vector  fields  u  and  v   with  the  aid  of  the 

■¥ 

vector  field  A, namely, 

u  =  IkA  -  rr-    grad  div  A 
(7.5)  .^^ 

V  =  curl  A, 

and  demonstrate  that  u  &nd  v  satisfy  the  equations  (7.2).  Indeed 

curl  u  =  ik  curl  X  =  ikv, 
i.e.,  the  second  equation  (7.2)  is  satisfied.  Furthermore 

curl  V  +  iku  =  curl  curl  X  -  k  A  -  gfad  div  X 

=  -4  a"-  k^  r=-ik|-  , 
i.e.,  the  first  equation  (7.2)  is  satisfied. 

We  thus  recognize  in  (7.5)  solutions  of  the  steady  state  equations 
(7.2).   Of  course,  we  cannot  be  sure  that  this  particular  soliition  introduced 
somewhat  arbitrarily,  represents  the  periodic  field  which  will  "be  established 
with  increasiiig  time  under  the  influence  of  the  enforced  currents  (7.1).  All 
we  have  shown  is  that  the  integral  fsatisfies  Maxwell's  equations  (7.2).  Actually 
however,  this  is  the  case  as  we  shall  verify  in  a  follovring  chapter. 

7.3.  Let  us  now  assume  that  enforced  currents  are  confined  to  a  small 
neighborhood  of  the  point  ^  =  ^  =  ^  =  0.  Let 

f-  1 ^  "m      for  t^  +  n^  +  i;2  <  /2 

=  0  for  ^^  +/,^  +^2  >  /2 

where  h   is  an  arbitrary  constant  vector  and  a   is  arbitrary  but  small. 
We  find  for  O  "-►  0  that 

->»     .^  J-^  ,/~p p p— 

(7.7)  A  =  -  M  -j^  ,  where  r  =Vx'^  +  y  +  z 

By  inserting  A  into  the  general  result    (7.5)  we  get  Hertz*   solution 
for  the  radiation  of  a  harmonic  dipole  of  moment     M  in  a  vacuum,   naiaely, 

V  =  -   (grad  rx  M)   (|  -     -^)  e^^ 

(7.8)  '^^ 

1 
u  =  -  rr  curl  V, 
ik 
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We  note  from  (7.8)  that  the  electromagnetic  field  of  dlpole  soiirce 
can  "be  written  in  the  following  form 

('•5'  Ikr      1 

where  the  vectors  Uj^   and  V^^   are  given  "by  the  relations 


(7.10) 


^0  ^  "  r  ^^^^'^  ^y''^^'  ^0  ^  "  ^^^'^  ^^^c 


T-  =  +  —  (grad  rXM),    ^i  ~  ~  S^ad  r  xV.  -  curl  7 
r 


TJ^  =  -  curl  Y^   , 

Since  these  vectors  are  independent  of  k  we  conclude  from  (7.9)  that 
with  increasing  values  of  k  =  "5"  •»  ^.®.i  ^°^  small  vmve  length  ^,  the  field  is 
asymptotically  given  hy  the  vector  field 

u=U  e^^^ 

(7.11) 

_   ikr 
V  =  V  e 
0 

We  also  note,  for  a  fixed  value  of  k  hut  large  values  of  r,  that 

the  main  contrihution  to  the  field  comes  from  the  expressions  C7.1l)» 

7,lV.  The  vectors  U  and  V  v/hich  determine  the  asymptotic  "behavior  of  the 
field  for  large  values  of  k  or  r  are  explicitly  given  hy  the  expressions 

U^  =  -  p  (grad  rXM)  X  grad  r 

(7.12)  1         ^ 

7    =  -  -   (grad  rXM)  . 

Because  €  and  M  are  constant  in  our  space  the  surf  aces  >!'' =  r  =  const, 
represent  a  system  of  wave  fronts.  We  may  therefore  interpret  the  unit  vectors 
grad  r  =  ^  as  wave  normals  and  write  (7.12)  in  the  form 

u^  =  -  p  (^xl^)  /.  t 
^7-13)  y  =-1(^kU)  . 
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Ve  conclude  that 


so  that  U  and  V  are  tangential  to  the  wave  fronts  and  normal  to  each  other. 

0        0 

Furthermore  from  (7.I3)  it  follows  that 
(7.15) 


pxV  +  U  =  0 
■^0    0 


PXU  -  V  =  0   , 
•^'^  0    0       • 

and,  since  p  =  grad  r  is  constant  along  a  straight  line  through  the  origin,  it  follows 

from  equations  (7.12)  that 

(7.16)  f   ° 

If  we  write  these  last   ecpiations   in  the  form 

dr  2  "  ^      0 

dV 

where'>K=  r  and  ^Y  »   accordingly,  is  —  we  recognize  "by  comparison  with  (5.1?) 

that  U  and  V  satisfy  the  same  relations  as  those  derived  for  electromagnetic 
0      0 

signals,  i.e.,  for  vectors  E*  and  H*  on  the  toundary  of  a  region  of  penetration. 

This  significant  fact,  that  the  amplitudes  of  periodic  waves  of  small 
wave  length  are  propagated  according  fco  the  same  laws  as  electromagnetic  signals, 
will  be  established  shortly  in  general* 

7.5.  We  consider  the  radiation  problem  for  a  time  periodic  current  distri- 
bution (7.1)  in.  the  case  of  a  non-homogeneous  medium.  The  steady  state  solution 
u,v  must  satisfy  the  differential  equations  (2.8),  namely, 

cu^l  V  +  ik6  u  -  -^u  =  -  ikl" 

c 

(7.18) 

curl  u  -  ikx  V  =  0. 
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We  assume  that  £,^,^  are  continuous  functions  of  x,y,z  and  that  they  have 

sectionally  continuous  derivatives. 

We  ohtain  the  radiation  from  a  dipole  in  a  similao"  manner  as  tefore 

"by  choosing  the  special  vector  field  given  by  (7.6)  as  our  source.  Since  this 

source  is  zero  except, at  x  =  y  =  z  =  0,  at  all  other  points  the  field  it  gives 

rise  to  must  satisfy  the  equations 

Utt 
curl  V  +  ik  £•  u  -  — cru  =  0 

(7.19) 

curl  u  -  ikyU  V  =  0. 

In  the  neighborhood  ofx=y=   z=0,   the  radiation  field  will  have  the  same 

singular  behavior  as  the  radiation  solution  for  a  homogeneous  space  in  which 

€  =    fe    .  /^  =    /^    1    <r=  0^     are  constant. 

We  consider  next  the  wave  fronts ■%(/*( x,y,z)  =  const  of  our  medium 

around  the  point  x  =  y  =   z  =  0,     We  know  that  the  wave  function  "/^satisfies  the 

differential  equation  (U.I5),   namely, 

(7.20)  ^l^^->^l     =€/^    =u2. 

In  view  of  the  results  for  the  homogeneous  medium  (7.9)  in  particular, 
we  can  make  the  attempt  to  represent  the  solution  u,v  of  Maxwell's  equations  (7.I9) 
by  expressions  of  the  fonn 

Let  us  derive  the  necessary  conditions  for  the  coefficients  U^  ,  V^  in  such  a 
development. 

We  introduce  the  infinite  series  (7.21)  into  the  differential  equa- 
tions (7.19)  BXii.  gToa-p   terms  in  accordance  with  the  power  of  k.  We  may  then 
equate  to  zero  the  coefficient  of  each  power  of  k  and  we  find  that  the  coeffici- 
ents U-^  and  7j^    must  satisfy  the  conditions 


(7.22) 

and 

(7.23) 


grad-i/'X  V     +  6U    =0 

'  o  0 

&TaAyfr}fi  U^  -  /A  7^  =  0 

grad>|^X  V^  +  €1^,  =  -  (curl  Y^_^-  ^<r  I{^_^) 

grad>f^;(  V^-  M\=  -  curl  ^_^ 


-Uo- 


7, 7. Consider  first  the  quantities  U  and  V  which,  in  view  of  (7, 21) 

determine  the  electromagnetic  field  for  small  wave  length.  \Je   have,  intro- 
ducing the  time  variation. 


(7.2U) 


and  we  recognize  that  wave  fronts  >^=  const,    represent   surfaces  at   equal  phase 
for  these  waves. 

The  equf.tions   (7.22)  are  fomally  identical  with  the  conditions   (U.12) 
for  discontinuities   on  wave  fronts.     By  substituting  one  equation  of   (7.22)  Into 
the  other  we  see  that   the  two  equations  can  he  satisfied  "by  vectors  U     f^  0  and 

V     ^  0  only     ifyyie  a  solution  of  the  equatiouAl/"     +'\P^     '*^^f^     =  n  ,     A  develop- 
ment of  the  form  (7.21)   thus  is  possible  only  if  the  function'y^satisfies  this 
differential  equation.      The  difference  between  the  equations    (7.22)  and  the  equa- 
tions   (U.I2)  for  wave  front  discontinuities  is  that  U     and  V     are  in  general  com- 

0  0 

plex  vectors. 

We  conclude  readily  from  the  equations   (7.22)  that 

(7.25)  U     •   greLd.yfr=  V       •   grad'/^=  U     •  V     =  0 

i.e.  U  and  V  are  tangential  to  the  wave  fronts  and  normal  to  each  other. 
— o 0 ■ 

We  consider  next  the  equations  (7.23)  in  case  l/=  1.  We  have 

grad-^AV,    +  €  U.    =  -   (curl  V     -  —  <rU   ) 
^  •  1  1  0       c  o 

(7.26)  , 

grad>f^XU^  -   >UV^  =  -  curl  U     . 

These  equations  are  formally  identical  with  the  equations  (5.7)   for  the  vectors 
E*   ,   -  —  H*   ,  E*,  &  H*  on  a  wave  front.     By  the   same  argument  we  conclude 

C        X  C        V 

that  they  are  solvable  only  if  the  right  sides  satisfy  certain  conditions. 
The  former  method  of  deriving  these  conditions  can  be  repeated  literally.  We 
obtain  the  result:   the  amplitude  vectors  U  and  V  satisfy  the  following  differ- 
ential equations  along  the  rays: 


which  are  identical  with  (5. 17).  Hence  in  the  case  of  the  dipole  in  the  non- 
homogeneous  medium  it  is  again  true  that  for  small  wave  lengths  the  "behavior 
of  the  steady  state  field  is  the  tehavior  of  signals  on  a  wave  front  vAiich  in 
turn  is  the  "behavior  of  geometrical  optics, 

7. 8. Our  treatment  of  the  steady  state  dipole  field  has  thus  far  assumed 
that  €  ,  /4  ,  <r  are  either  constant  or  continuous.  The  asymptotic  representa- 
tion assumed  for  the  steady  state  field  is  relatively  simple  in  these  cases,  the 
more  complex  "being  (7. 21)  for  the  continuous  case.   If  € ,  JU  t        ^^^  ^   are  dis- 
continuous then  the  asymptotic  expansion  "becomes  more  complex. 

Discontinuities  in  6  ,  /^  ,  and  C  are  generally  caused  "by  the  presence 
of  refracting  surfaces.  We  therefore  consider  the  case  that  our  medium  conteine 
a  nam"ber  of  refracting  surfaces  X^  on  which  the  functions  6,yW  ,  CT  are  discontinu- 
ous. We  assume  that  a  point  source  is  located  at  the  origin  so  that  the  vector  field 
g(x,y,z)  is  zero  except  in  the  neigh'borhocd  of  the  point  x  =  y  =  z  =  0,  Let  us  also 
assume  that  the  refracting  surfaces  are  either  closed  or  extend  towards  infinity 
and  that  the  space  is  divided  into  a  finite  num"ber  of  domains  in  which  €  ,/4  ,cr 
are  continuous.   On  account  of  multiple  reflections  we  o"btain,  in  each  of  these 
domains,  infinitely  many  sets  of  wave  fronts "^^1^=  const.   In  view  of  the  former 
results  we  attenrpt  to  represent  the  steady  state  solution  u,v  in  form  of  an  infin- 
ite series 

u=      ru,e"^>^ 
(7.28) 

V  =     Z_  Va(  e      '^ 

ik"^/^  ik'v/c 

where  each  individual  term  U^t    e  and  Vo<.  e       '^  satisfies  Maxwell's   equations. 

The  representation  (7.28)  is  a  generalization  of  the   representation  (J. 21)  in  that 

it  allows   for  many  sets  of  v/ave  fronts  possi"bly  an  infinite  iKim"ber,  passing  any 

given  point    (x,y,z)   in  space  at  which  the  steady  state  field  u,v  is  "being  computed. 

Each  Uv     and  "7^     is  assumed  to  "be  an  infinite  series  of  the  form  given  on  the  right 

side  of  (7.21),   that   is,    each  l^     (and  Vot  )  is  of  the  form    ^^  r^^vv 

To  determine  the  various  Uot  ,v  and  It*,  tv  it  is  necessary  to  use  the  usual 
electromsignetic  "boujidRry  conditions   at  the  surfaces  X     ,      It   is  not  difficult   to 
determine  the  field  emanating  directly  from  the  source  and  transmitted  through  any 
one  of  the  surfaces  X    .      In  principle  it  is  also  pos8i"ble  to  construct  the   total 


small  wave  length  field,    i.e.,    to  find  the  various   fields  which  correspond  to 
the  infinitely  many  internal  reflections.     Fortunately  an  e:act  knowledge  of 
these  secondary  fields  is  not   often  of  practical   interest.      In  optical  instru- 
ments  one  tries  to  reduce  these     reflections  to  a  miniimiin  absorption  since  they 
disturh  by   "ghost  images"   the   image  formation  of  the   instrument. 

Further  discussion  of  the  treatment   of  the  full  steady^tate  field 
arising  from  a  dipole   source  will  not  be  given  here,    for  the  purpose  of  this 
article  is  merely  to  illustrate  the  general  theory  of  asyinptotic  development 
of  steady  state  fields,   which  will  be  begun  in  the  next  article, 

8.     Application  of  Laplace  Treuasform  Theory  to  Itezwell's  Equations. 

8.1. This  article  will  present   some  general   theory  on  the  solution  of 
Maxwell's   equations  when  the  radiation  field  is   due   to  an  arbitrary  source.      It 
is   this   theory  which  will  be  applied  in  the  next  article  to  deduce  the  asymptotic 
expansion  of  the  stead;/'  state  field  and  through  this   the  relation  of  geometrical 
optics   to  the  steady'  state  field  for  any  given  wave-length. 

The  essence  of  this  article  is  one  basic  theorem.      In  effect  the 
theorem  accomplishes  the  following.      The  problem  of  solving  Maxwell's  equations 
for  E(x,y,z,t)  and  H(x,y,z,t)   is  attacked  by  introducing  the  Laplace  transforms 
with  respect  to  time  of  E  and  H.      These  transforms,   W  and»fl»which  are  independent 
of  time  are  shown  to  satisfy  the  simpler  fonn  of  Maxwell's  equations  which  one 
obtains  uhen  the  source  is     harmonic.      It   is  then  possible  to  define  two  vector 
functions  A  and  B  in  terms  of  W  andtTL  ,   from  \Ai.ich  B  and  H  can  be  obtained  by 
simple  differentiation. 

What  the  theorem  accomplishes  is  that   it  furnishes  the  expressions 
(8.10)   for  E  and  H  in  terms  of  integrals   in  which  the  variable  t  enters  as  e        , 
Hence  the  expressions  hold  regardless  of  the  character  of  the  source,   that  is, 
whether  it  be  discontinuous  or  continuous.     Of  course  implicitly  the  character 
of  the   source  nnist  be  involved  in  the  integrals   through  the  Laplace  transformation. 
The  usefulness  of   the  theorem  will  be  apparent   in  the  next  article. 

8, 2, We  alter  the  problem  of   solving  Maxwell's  equations    (2,2)  by  introducing 
the  Laplace  adjoints  ^^ 

W(x,y,z,r  )  =   J      E  e"  ***  dt 


(8.1)  ^(x,y,z,r)  =  J     H  e'^* 


f  (x,y,z,r)=    I      F  e  '^ 


r 


dt 


dt. 
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It  is  a  fact   of  Laplace  transform  theory  that  at  any  fixed  point 
x.y.z,   where  E,H,    and  F  are  bounded  and  sectionally  smooth  in  t   then  W.^TL, 
and  J    will  be  analytic  functions  of   the  complex  parameter  /* =  of +  iyg   which 
are  regular  in  the  half  plane  ot  >  0. 

Let  us  now  assume  tliat,    at  a  fixed  point  x,y,z,    one   or  more  of  the 
vectors  E,  H,   or  F  are   discontinuous   at   the  times 

Any  one  function  t^  represents  the  functional  dependence  of  time  on  position  in 
space  so  that  as  one  passes  from  one  point  to  another  t^^    will  vary  gradually. 
At  any  fixed  point   of  space  there  may  be  many  functions  -^^  which  represent  differ- 
ent discontinuities,    resulting  perhaps  from  many  separate  wave  fronts  passing  that 
point  at  different   times  tex    ,    the  wave-fronts  themselves  resulting  from  reflec- 
tions from  different  obstacles   in  the  medium.      It   is  also  possible  to  regard  the 
many  functions  "^^  as  branches  of  some  implicit  function,    ^(x,y,z,t)  =  0. 

The  hypersurfaces  ^f^  (x,y,z)  =  ct  divide  the  x,y,z,t   space  into 
separate  regions  in  which  E,  H,    or  F  are  sectionally   smooth.      From  the  discontinu- 
ity conditions  (3.9)  which  apply  to  any  hypersurfece  ^(x,y,z,t)  -  0  it  follows 
that  E,H,   and  F  must  satisfy  on    W)  ^^lr{r.,y,z)  -  ct  =  0,   the  conditions 

grad^X  [h]  +    [ce]     =  -  [f] 
(8.2) 

grad-vfc  X  [e]  -    ^^hI     =     0. 

If  in  particular  £  ,  U  and  F  axe   continuous  on  such  a  hypersurface  "^fc  =   ct  then 
it  follows  by  the  argument  used  to  obtain  (U.I5)  that  "^  (x,y,z)  must  be  a  solu- 

P  •   P  P 

tion  of  the  eiconal  equations  ^1^     "*"  ^y  "'"'^    =  ^/^    .    i.e.,  VS  =  ct  mast  be  a  set 
of  wave  fronts.      In  addition  to-  such  wave  fronts  other  hyper  surfaces  may  occur 
v/hich  are  caused  by  the   discontinaities  of  €  ,//  or  of  the  vector  field  F. 

We  now  consider  a  fixed  point  x,y,z  at  which   ^  ,/M  and   O"  are  continu- 
ous,  that  is,  we  exclude,    for  exajnple,  points  on  a  refracting  surface.     By  apply- 
ing integration  by  parts  to  relations   (8.1)  we  obtain". 

erw  =   ZT  USJ^     e"  ^**+e  f\   e-^*   dt 
^  o 


(8.3)  ><''■«•=    1^    fy«<     ^"'^'"'•/m/  H,  e-*"'  dt 

ri^  :^  M,  e-  »•'■ .  Sj.  -- "'  ^'  • 

9.  To  obtain  this  result  we  must  regard  the  integral  from  0  to  00  in  (8.1)  as  the 
sum  of  the  integrals  from  0  to  t^  ,  from  ti   to  tg   ,  etc.,  and  lastly,  the  inte- 
gral from  ta.   to  00.  The  toe  are   the  values  of  t  at  which  E,H  and  F  are  dis- 
continuous at  a  given  point  (x.y.z).  Alternatively,  v/e  could  introduce  the 
Stieltjes  integral  in  (8.3)  and  this  result  would  follow  at  once. 


-HV 


We  now  return  to   (8.1)  and  regard  it   as  the  sum  indicated  in  footnote  9.     By 
differentiating  under  the  integral  sign  in  (8.1)   and  rememhering  that  the  dis- 
continuities  in  E  and  H  are  normal  to   ^Q  (x,y,z)  we  get. 


(8.U) 


curl  V  =  -  i-Z^grad-^A  /[E^e"  ^W  ^  curl  E  e"  ^*  dt. 
curl  A  =  -  i-ZT  grad^CX  CnLe"  '^^fj^rl  H  e"  ''*  dt. 


Thus   it  follows  "by  Maxwell's  equations   (2.2)  and  the  disco  ntlmiity  relet  ions 
(8.2)  that  W  and»fL  satisfy  the  equations 


(8.5) 


curlA-     -Cfe  v-lj!(rW=>y$ 
c  C  *" 

curl  W  +     _llyU^=  o. 


The  Laplace  adjoints  V   and ^ of  S  and  H  thus  nust  he  solutions  of  these  equations 
which,  at  any  fixed  point  x,y  z,  are  analytic  functions  of  the  com-olex  t»araineter 
<*  and  are  regular  v/ith  half  plane  Ry]  >  0. 

The  continuity  of  the  vectors  grad  ^  X   E  and  grad^  X  H  on  any  surface 
(^(x,y,z)  =  0  of  refraction  yields  the  corresponding  continuity  conditions 
(8.6)  grad<|)>tlw]  =  grad0xlA)=  o 

on  any  surface.  This  follows  "by  crossing  "both  sides  of  (8.1)  "by  grad  ^  and  putting 
grad 0  under  the  integral  sign  on  the  right  side. 

Under  sufficiently  restrictive  conditions  on  E  and  H  it  is  possible  to 
invert  the  Laplace  transforms  "by  means  of  the  usual  inversion  theorem,  namely, 

J  W  e  ^*  dr 


(8.7) 


E  = 


H  = 


F  = 


1 

2ni 


1 
2ni 

1 
2ni 


Jj^e^* 


dr 


I  §e^'   6.r 


-»-o( 


The  path  of  integration,  L,   can  "be  aay 
straight  lineO(  =  c/     >  0  T3arallel  to 

0 

the  imaginary  axis.     However,   we  shall 
not  impose  such  conditions  hut   shall 
deal  directly  with  the   functions  W,*^ 
and  $    . 


Pig.  13 
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S,3»   ^^6  problem  of  constructing  solutions  of  the  original   radiation 
proljlem  (2,2)  can  be  replaced  "by  the  problem  of  finding  solutions  of  the  equa- 
tions  (S,5)  which  satisfy  the  boundary  conditions   (8.6)  and  regularity  condi- 
tions in  the  half  plane  H  \_'''*)>  0.      Indeed  we  have  the  following  Theorem. 

1)  Let  F(x,y, z,t)  be  a  sectionally  smooth  and  uniformly  bounded  vector 
field;      i.e.,  |F|  <  C   for  all  x.y.z.t;      also  let  ^  as  defined  by  (8.l)  be  its 
Laplace  adjoint, 

2)  Let  V  and^be  solutions   of  the  differential  equations   (S.5)  which  satis- 
fy the  boundary  conditions   (8.6)  and  the  condition  that  they  be  continuous  and 
have  continuous  derivatives  in  any  closed  subdomain  D  of  the  x,y,z  space  where 

^  ,  JU  ,   and  (Tare  continuous. 

3)  It  is  further  assumed  that  to  any  such  subdomain  D  of  the  x,y,z  space 
there  exists  a  positive  function  M(P  )  such  that  J  H(/0)  d^  exists  and  such 
that  ° 

in  D  and  for  all  f=0/+  lytf  in  the  half-plane     «^   />  0. 
k)  Let  the  integrals  A  and"^,    defined  by 

(S.9) 


^-     arii     -?2    ^        ^^ 


'L  r2  " 


represent   two  vector  fields  which  have  continuous  first   derivatives  and  sec- 
tionally continuous   second  and  third  darivatiwes    in  D  and  for  t  ^  0, 

It  follows   that  the  vector  fields 

(g.lO)  E  =     ^-4    and    H     =     ^  — 


9t2  dt2 

satisfy  Maxwell's  equations   (2.2),   the  discontinuity  conditions   (8.2),   the  con- 
ditions,  grad<^x[lj=  grad<f>^[E7=  0  on  the  refracting  surface  <f){x,y,z)  =  0, 
and  E  and  E  are  zero  for  t  <  0. 


-U6- 


Proof.      Consider  a  closed  domain  D  of  the  x,y,z  space  in  which  €.,  A*  ,  <r  are 
continuous  and  assume  that  a  part   of  its  "botindary  surface  coincides  with  a  sur- 
face of  refraction.      Our  assumption  ahout  M  andJL  means  that  "both  vector  fields 

attain  definite  limits  if  a  point  P 
of  1  is  approached.  Since  the  inte- 
grals (8.9)  are  uniformly  convergent 
in  D  and  for  t  2  0  by  the  hypothesis 
of  the  theorem,  it  follows  that  also 
A  and  B  attain  definite  limits  given 
hy  the   integrals   (8.9).     As   in  the 

proof  of  (8.6)  it  now  follows  that 

^^S.    1^  the  condition  grad^J^[1tJ=  grad0X^j=O 

then  is  a  direct  consequence  of  the   corresponding  conditions   for  V  and/t  on  T 
The  existence  of  the  second  time  derivatives  of  A  and  B  and  the  interchange  of 
order  of  differentiation  in  (S.IO)  leads  to  the  conclusion  that  grad^*^]  = 
grad4>X  Lh]      =  0. 

We  next  conclude  from  (S.8)  and  (S.9)   that 


(S.ll) 


bI^ 


Md 


? 


,^1   .       6 


"  0 

where  O^    ^  o  >  0  is  the  arbitrary  abscissa  of  the  path  of  integration,  L.      If 

t   <  0  then   (A/ and  {%\  are  smaller  than  an  arbitrary  small  number  i.e.     A  =  B  =  0. 

Hence  by  (8.10)  E  =  H  =  0  for  t  <  0. 

Tor  the  proof  that  the  Maxv/ell  differential  equations    (2.2)  are  satis- 
fied note  that  the  hypothesis 

IfI  <  C 

inplies  for  the  Laplace  adjoint  ©  that 

*'o 

i.e.,  W  is  uniformly  bounded  in  half  plane    O^^  /  >  0.     It  also   follows  from 
(8.7)  and  a  basic  theorem  on  the  Laplace  transformation  that    the  Laplace  adjoints 


of  the  iterated. integrals 

■t 

F(T)  dr     and 
o 


/ 
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are  given  "by  the  expressions     _*.    and  -^     and  hence  by  Laplace's   inversion 
taeorem  that 

L 
Replacing  the  left   side   of  this  second  equation  hy  an  equivalent   expression  gives; 

<^-^'>  /;u-r).(r)dT=  ^/-£/«.r. 

We  now  consider  the  vector  fields  a!.    and"!!    defined  hy 

T    =       1/   J^e^^d/- 


(S.lU) 

L  r3 


Since  by  formal  differential  with  respect  to  t  unifonaly  convergen-^  integrals 
are  obtained  we  conclude  that     -^ 

A  = 


a  t 

(8.15)  »-  9^1 

^  ~  ©t 

Moreover  if  in  (S.5)  we  divide  by  -f  .  multiply  by  e   ,  and  integrate  we  obtain 


(8. 16) 


J        Y" 3  °  •'J^^  '^       •'^3  ° •'  r^ 

fcurlJLe'^dr=-^f^e^*dr      , 

i.e.,    the  curl  operation  under  the  integral  sign  yields  uniformly  convergent 
integrals.     Hence  by  (8.IU) 

T  -*►    _       1      fcurl  W       n     ^  - 

eurlt    =     ^^f^Hl:^     e'*     dr 
1        2ni  Jj^   1%  3 

and  by  (S.15),    (8.I6)  and  (S.!"?)  weJiave 

-     eurl  I;   -     ft    |i    J^f    =  I  J  '  (t-  r)F(r )  dr 
1        cgt        *^-i-      ^  ''o 

(S.17)    ■  *     >f  a% 

curl  Ai    ■►  -^    -^-i-    =     0. 
^         cat 
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In  any  domain  of  the  x,y,z,t  space  where  A  and  B  have  continuous  time  derivatives 
to  the  third  order  we  may  differentiate  (8,17)  three  times  with  regard  to  t.  We 

ottain  immediately  the  desired  differential  equations  (2.?)  for  E  =  -      and 

9^  a  t^ 

In  order  to  verify  the  discontinuity  conditions  (8.2)  let  us  differen- 
tiate (8.17)  twice  with  respect  to  t.   It  follows  ^hat 


(S.18) 


eurl   £|  -  £  E-  ^(T^  =  i 
9 1      C       C    9  t      c 

curl   *^  -*•  ^  E  =  0. 

at         c 

.2^  -2|- 


3  A  S 

Let  us  now  assume  that  E  =  ^  and  H  =  ^ — ^  are  discontinuous  on  a  hyper  sur- 

^*  at^  ^         -. 

face>I^(x,y,z)  =  ct,  which  lies  in  the  domain  D  for  t  >  0.  Since  ^4  and  SJ.  are 

9  t      ©  t 

continuous  in  D  "by  hypothesis,  and  the  equations  (8.18)  are  satisfied  on  both  sides 
of  the  surface  it  follov;s  from  (8.18)  that 

[-1  If]  -  f  W  =  I UJ 

(8.19)  ^ 

[--1  It  J  *  f  W  =  °. 

where  curl  B.  means  curl  B. (x,y,z,t).   Let  B*  he  the  value  of  3  on  the  surface,  that 
is,  B*  =  B(x,y,z,  -^  (x,y,z)j  .  Let  ^-^  =  B*  mean  — r  and  suhstitute  •^^(x,y,z) 
for  t  after  differentiation.  Then,  for  example, 

S      ■*-  »•*■■♦"  ^y 

ay    t       ay    t    *tt     c 

Prom  this  it   follows   that 


(8.20) 

curl 


C 

the  two  equalities  hold  only  on  the  surface>^=  ct.  The  left  sides  are  zeio  since 


where  -'—   is  substituted  for  t  on  the  right  side  after  differentiation  and  hence 
c 
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ISA*')       f3B*7 
St~J  ~    ^aT~^  ~  *^   '      ^^'^'^^^  ^y  ^^^   definitions  of  E  and  H  and  "by  (8.I9)  and 

(s.;?o). 

grad^A^H]  +  €[bJ  =  -        [f] 

(S.21) 

grad-ZoilE]  -  /<fH]  =  0 

and  these  are  the  conditions  (8.2)  if  we  represent  the  hyper  surface  in  the  form 
(u=    j^-   ct  =  0.   This  completes  the  proof  of  the  theorem. 

9.  The  Asymptotic  Development  of  the  Steady  State  Field. 

9,1. This  article  will  employ  the  theorem  Just  estahlished  to  arrive  at 
an  asymptotic  expansion  of  the  steady  state  field  created  by  harmonically  oscilla- 
ting charges.  The  expansion  to  be  derived  is  a  series  in  which  the  variable  is  the 

wave-length  K  ov  —    where  w  =  2n  tiHies  the  frequency.   The  coefficients  of  the  powers 
of  X  have  the  following  meaning.   One  defines  first  the  pulse  solution  of  HaJtwell's 
equations.   This  is  the  field  whi^ih  is  set  up  by  a  source  which  is  0  for  t  <  0  and 
which  for  t  ^  0  is  1  in  amplitude;   that  is,  at  t  =  0  a  charge  having  the  same 
spatial  distribution  as  the  harmonically  oscillating  charge  is  suddenly  set  up. 
The  resulting  field  satisfying  Maxwell's  equation  ie  the  pulse  field  for  this 
charge  end  is  denoted  by  E  and  H  ,  At  any  point  in  space  this  field  has  dis- 
continuities with  respect  to  the  time  variable,  the  discontinuitj  s  arising  either 
from  the  distribution  of  the  sovirce,  from  the  passage  by  the  point  of  the  v;ave-front 
emanating  from  the  source,  or  from  discontinuities  in  €  ,d4.   ,  and  (T  which  cause  re- 
flected wave  fronts  to  pass  the  point. 

The  essential  parts  of  coefficients  of  the  terms  in  the  asymptotic  ex- 
pansion for  u,  the  spatial  component  of  the  steady  state  E,  are  respectively,  the 
sum  of  the  discontinuities  of  E  ,  the  sum  of  the  discontinuities  in  the  first 
derivative  of  E  ,  the  sum  of  the  discontinuities  in  the  second  derivative  of  E  ,etc. 
The  same  s-oatements  apply  to  v,  the  spatial  factor  of  H,  and  H  . 

The  significance  of  this  result  which  is  embodied  in  formula  (9. 18) 
will  be  coamented  upon  at  greater  length  at  the  end  of  the  article.  It  will  be  seen 
shortly  that  the  theory  of  this  article  generalizes  the  result  noted  for  the  dipole, 
namely,  that  for  small  wave-lengths  the  amplitudes  of  the  steady  state  electromag- 
netic field  obey  the  same  mathematical  laws  as  do  the  amplitudes  of  the  vectors  on 
a  progressing  wave- front. 
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9.2  Let  us   consider  vector  fields  F  of  the  special  form 
(9.1)  ?=T(x.y,z)  f(t) 

where  "fCx.y.z)  is  a  given  vector  field  and  f(t)  a  giren  scalar  function.     Ve  do 
not  exclude  the  case  that  the  components  of  g  and  the  function  f(t)  are  complex 
functions   of  their  argument.     We  assume  however  that  g  is  sectional ly  smooth  and 
uniformly  "bounded  in  the  x,y,z   space.      Similarly  f(t)  shall  be  a  sectionally 
smooth  and  uniformly  hounded  function  of  t.     Moreover,    it  shall   satisfy  the  condition 

(9.=)  rj;f(oe-^'dti  =(<<■) 

in  any  half  plane  B.{f(^  (f  >  0, 

By  equation  (8.5)  of  the  preceding  article  the  Laplace  adjoints  of 
the  solutions  E  and  H  of  the  radiation  prohlem  of  article  2  satisfy  the  differ- 
ential equations. 

curl/I-     -^eW-     ^O-W     =     i?  rf  ^  f(t)  e''"*  dt 
c  c  c  J^ 

^^'^^  curl  W  *    ^rSl=  0. 

c 

9.3  Consider  in  poxtlcular  the  case 


(9.^)  f(t)  =  9(t) 


fO,  t  <  0 
1  t  >  0. 


Tliis  means  physically   (see  article  2.1)  that,   at  the  time  t  =  0,    suddenly  a  charge 
distrihution  /0=  -  div  g  is  created  in  the  space.     This  act  generates  an  electro- 
magnetic f ield  E  (x,y,z,t),  H  (x,y,z,t)  which  with  t  -►   00  mast  approach  stationary 

values,   namely  the  stationary  electromagnetic  field  of  a  charge  distrihution 

In  particular  H  (x.y.z.t)"^  0,     Let  us  call  this  narticular  solution  E   ,  H     the 
■  0  *  00 

pulse  solution  of  the  medium  and  of  the  charge  distrlTjiition^Z-  div    "^  . 

For  this  f(t),  the  right  side  of  the  first  equation  (9.3)  is  -  . 
Let  us  now  assume  that  the  associated  equations  for  the  Laplace  adjoints  W  and 
f\       of  E     and  H  ,   respectively,  namely, 

0  c  0  C  0  c 

(9.5)  yv      ^ 

curl  W       +    -i-iWXi.     =  0 

0  c   '       0 
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have  a  solution*/!,  ,  W     which  satisfies  the  premises  of  the  general  theorem  of 
article  8.      Then  by  the  theorem  of  that  article    (see  fonmilas  S.IO)  B     and  H 

0  0 


are  given  "by  the  integrals 


E    =    '•' -=-    I      — ^    e"^  dr 

\       ^     2     2ni  \     y2  '^^ 

(9.6)  ^'  ^    ^ 

a  ^     1   r  »A.o    rt  .^ 

\     ^  2Fi  J,  72  «    ^^ 

That  this  S     and  H     are   identical  with  the  pulse   solution  follows  from  a  imique- 
ness  theorem  on  the  solutions  of  Maxwell's  equations. 

Also,    since  W^  is  the  Laplace  adjoint  of  S     and««4«  ,   the  adjoint  of 
o  00 

H   ,   it  follows  by  a  basic  theorem  of  Laplace  transform  theory^*^  that 


L   I  0 

L   »  "0 


"-        A= 


This  result  tells  us  what  the  inverse  L^lace  transforms  of  — %  and  — -  are. 

yS      y2 

9,1+  Consider  next  a  more  general  function  f(t)  which  satisfies  the  condi- 
tions (9.2).  In  view  of  (9.5)  a  solution  ¥,A  of  the  associated  equations  (9.3)  is 
immediately  given  by 

■'  =  ^^r.   I         f(t)  e"  ^*  dt 
o   ./ 

(9.S)  ° 

jci.  ^rsuf"^  ^^*)  e-''*  dt  . 
"0 

On  account  of  (9.2)   this   solution  satisfies  the  premises  of   the   theorem  of  the  pre- 
ceding article.     By   (8.IO)  and  (9.8)  the  electromagnetic  field  E,  H  which  belongs  to 
F  =  g  f(t)  thus  is  given  by  the  Integrals 

10.  The  theorem  applied  is  the   "Paltung"  theorem.      It  can  be  stated  thus:      If 

y    Fj^(t)e"^  dt  =  ft  <r  )  and  /     T^it)e'''^^  dt  =(i^{r  )  then  -^f  ^i(X)(}^{r)^^  dr  = 
*  ■  L 

J*  I[(y)  ^(t-y;^.v.      In  our  sp-plication  5|.(t)  =  Bq  and  ^2^^)  =  *.   T^^en  Gi(Y")=Wo  and 
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The  uniform  convergence  of  these  integrals  maker,  it  possible  to  rewrite  these  inte- 
grals as  follows: 


(9.10) 


It   is    readily  seen  that    the  new  nniltiple   integrals   converge  uniformly 

in  If  and  t   and  thus    the   order  of   intet;ration  can  "be   interchanged;    that   is,   we   shall 

integrate  \/ith  respect   to  y*  first.      If  we  rewrite  7*  as  t  and  again  apply  the  Faltung 

theorem  to  the  functions   f(t)  and  J       (t-  1*)  ^  (T*  )    d.T   we  obtain  through  use  of 

,  •'o  ° 

fonrala   (9.7). 


(9.11) 


9 1-  -0 


t-T 


H=    ^-^J    f(r)/       (t-r-r.)  H  (r  )dr,  dr 

^  t^     o  0  loll 

and  finally  by  carrying  out  two  cf  the  differentiations,  treating  t  as  the  para- 
meter and  applying  the  theorem  for  the  differentiation  of  an  integral  with  respect 
to  a  pn.raineter,  '-'c  get,  ^.^ 

''  It  J 

(9.12)  ^t 


-.  .    f(r)  E  (t-r)  dr 

S^   0  ° 


H  =  -|^4  ^^'^'^  H^(t-r)  dr  . 

This  simple  result  allow?  us  to  determine  the  solution  E,  H  for  r.ny  source  of  the 
type  F  =  g(x,y,z)  f(t)  if  the  pulse  solution  for  the  case  F  =  g(y.,y,z)  'Q  (t)  is 
knovm.   The  formulae  (9.12)  are  'oiown  as  Duhamel's  Principle  in  Heayi^ide's  opera- 
tional calculus, 

9.5  IXihanel's  principle  establishes  in  particular  a  relationship  between 

the  pulse  solution  Eq,  H  and  the  solution  of  the  steady  state.   In  fact  consider  a 

harmonic  oscillator 

*  ■*,  X  -i<*'t 

(9.13)  F  =  elx,y,c)  e 

so  that  f(t)  =  e~     .   3y  making  the  change  of  variable  7*  for  t  -  T*  in  formal «i 

— ''cot 

(9.12)  and  using  f(t)  =  e  "^         we   find  tlriat  the   electromagnetic   field  associated 
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with  this   oscillator  is 


«9  L  "^  0 


^-  at  -        j^  -o- 


H=     It     ^-"^I     H„(r)e'"''dr 


or,  after  a  slight  trans format Ion 


E  =  B  (t)  -  E  (  oo)  +  e~^*^*r  E  (  oo)  -  ia)/    (a  -E(oo))e^**'  d.r] 

0         0  *■   O  J       0     0 

•0 

(9.1U) 

H  =  H^(t)  -  H^(oo)  +  e"^'***[H^(oo)  -  icof  {E^-E^{oo))  e^"^^    dth 

This  result  makes  it  evident  that,  with  increasing  time,  E  and  H  assume  the 


-iojt   „     -ion 
form  E=ue      ,H=ve      where 

Ti(x,y,s)  =  B^(oo)  -  iCOj   (B^  -  B^(oo);e    dT 
(9.15) 


(x,y,z)  =  H^(oo)  -  ia»r   (H^  -  H^(oo))  e^***''dr  . 

■'  n 


0 


In  these  fornrolas  E     and  H     axe  f-unctions  of  the  suppressed  varia"bles  z,y,   and  z, 

9,6  The  formulae  (9.15)  estahlish  a  simple   relationship  hetween  the  ampli- 
tudes u,v  of  a  time  periodic  electromagnetic  field  and  the  field  Eq,  H     of  a  sudden 
sigxial.      We  use  o\ir  result  to  determine  the  asymptotic  tehavior  of  u  and  v  for  large 
values   of  00,     We  integrate  hy  parts  and  remember   tliat  E     and  H     are  discontinuous 
at  the  hyper  surfaces  t  =  —  yCj  (3c,y,z),     Because  of  the  discontinuities  the  quantity 

E^e^***'*'  \     equals-23  [EqI,^  e^^'^^-yojllhe  full  result 


Is 


(9.16) 


u  =  B^(..y^z.o)  *  I  [.X  «''^    *f"  V^'"^^'' 


«         --   -  -'o 
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How,   from  (3,11 )  it  follows  that 

E^(x,y,z.O)=     -i    P(0)  =  -     I    g'U.y^z) 
H^(x,y,z,0)  =  0. 


This  yields 

the 

result: 

(9.17) 

u  » 

1 

o< 

9\ 

^ik^ 

H 
o 

< 

V 
V 

icar 

^ia;r 

dr 

dr 

If  also  the  derivatives     -  .        and    ^  are  sectionally  smooth-we  nay  again 

integrate  "by  parts  and  ohtain  the  representation 

u  =     -    i    1  ^        [2  L  e^^  -  J^rf  ^     e^^>^-  J- r^  e^««^dr 

romula  (9. 18)  is  of  course  a  complex  vector  equation.  Moreover  it  is  under- 
stood that  the  left  and  right  sides  are  functions  of  x,y,  and  z.  The  formula 
shows  that  the  main  contrihution  to  u,v  for  large  values  of  A>is  given  ly  the 
expression 

(9.19)  e   »  oT» 


<y  *•  O^ 

It  will  "be  noted  that  fomrilas  (9. IS)  give  an  exact  representation  for  u  and  v. 

However,  if  the  successive  higher  derivatives  of  i!  and  H  are  sectionally  smooth 

00 

the  integration  hy  parts  can  he  continued  and  the  successive  terms  of  the  asymp- 
totic representation  ohtained.  This  latter  view  of  the  result  will  prove  more 
useful  in  application. 

At  any  point  in  a  domain  where  ^,/d   ,<7"  are  continuous  and  where 
there  are  no  sources  we  have  for  small  \   that 

rik  V^ 

(9.20) 
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As  remarked  in  article  7  in.  cormection.  with  foramlas  (7.22)   to   (7,2U)  the 

ftmctions     rfiC^.y.z),    in  this  case,  must  te  solutions  of  the  equation 

2  2  2 

^^    *  yUy  "*"  ^/^      ~^M    .i.e.,  arast  he  wave  fronts.     The   factors  Uq^    ,   7^    are 

given  hy  the  discontinuities  Uot    =     [^Qitt    i   ^     "  [H  Jo^of  the  pulse  solution 
Eq,  H     on  these  wave- fronts.      These  discontinuities  thus  determine  directly   the 
asymptotic  hehavior  of  the  steady   state  solution  in  the  case  of   small  wave  length. 

In  general  we  can  use   the  formulae   (9.18)  to  construct  steady  state 
solutions  from  the  particular  solution  which  represents  the  field  of  a  sudden 
introduction  of  electric  charge  into  the  space.      The  formulae  make  it   evident 
that  the  higher  terms  of  an  asyr^totic  development   of  u,v  in  case  of  large 
values  of  k  are  given  hy  the   discontinuities   of  the   derivatives     ^ — ?.      SI — Q. 

^  ^    dt^ 

etc.  of  the  pulse  solution. 

The  "basic  result  (9.18)  is  useful  hoth  theoretically  and  practically. 
On  the  theoretical  side  it  says  that  if  we  have  full  knowledge  of  the  pulse  solu- 
tion of  Maxwell's  equations  due  to  a  particular  source,  or  if  we  at  least  know 
the  discontinuities  of  that  pulse  solution  and  its  derivatives,  \/e  can  express 
the  steady  state  due  to  a  harmonic  oscillation  of  that  same  source  in  terms  of 
the  discontinuities  of  the  pulse  solution  and  its  derivatives.   While  the  analo- 
gous principle  is  knovm  and  has  been  used  in  circuit  theory,  it  has  special  sig- 
nifitfance  for  electromagnetic  theory.  As  the  text  points  out,  at  a  point  x,y,z 
in  a  domain  in  which  6,  JU,   and  ^ are  continuous  and  in  which  there  are  no  sources, 
the  discontinuities  must  arise  from  wave  fronts  passing  this  point.  Of  these 
wave-fronts  one  would  he  the  wave  front  arriving  directly  from  the  source,  the 
transmitted  wave-front.  The  other  wave-fronts  will  arise  by  reflection  from 
discontinuities  in  €  , /^  ,   and  <r  ,  existing  at  other  points  in  space,  and  due, 
for  example  to  the  presence  of  obstacles.  Now  it  was  shown  in  earlier  articles 
that  the  behavior  of  E  and  H  on  these  wave-fronts  is  the  geometrical  behavior 
usually  furnished  by  geometrical  optics.  Hence  the  discontinuities  of  E-  and 
Hq  on  the  wave  fronts  are  precisely  what  geometrical  optics  supplies  as  to  the 
field  at  the  point  (x,y,z).  In  other  words,  the  first  terms  in  u  and  v  of 
formula  (9.18),  except  for  the  source,  is  essentially  the  geometrical  optics 
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description  of  the  steadi'  state  field.  The  theory  shows  then  that  for  very 
small  wave  lengths,  the  geometrical  optics  solution  is  a  good  approximation 
and  at  least  a  first  approximation  to  the  steady  state  field. 

The  result  (9.18)  while  useful  in  a  variety  of  applications  is  one 
fonmilation  of  the  asymptotic  expansion  of  the  steady  state  field.  A  new  formu- 
lation and  its  advantages  will  he  given  in  the  next  surtlcle. 

On  the  practical  side  the  asymptotic  expansions  (9.18)  can  he  used 
to  obtain  approximate  values  for  the  steady  state  field.   If  only  the  geometrical 
optics  representation  of  the  field  is  known  then  the  theory  tells  us  at  least 
the  true  character  of  the  approximation  in  a  much  ftiller  and  more  explicit  manner 
than  hitherto  known.   If  in  addition  the  pulse  solution  is  known  several  terms  of 
the  asymptotic  expansion  can  he  computed.  It  will  he  shown  in  a  later  paper  that 
the  pulse  solution  can  he  obtained  in  some  diffraction  problems.   Of  course  a 
fixed  nomber  of  terms  of  the  asymptotic  expansion  will  be  a  better  approximation 
for  smaller  wave  length. 

10,  Complex  Integral  Representation  of  the  Asymptotic  Development. 

10.1  The  asymptotic  representation  of  the  steady  state  electromagnetic 
field  furnished  by  formula  (9.18)  is  adequate  when  the  derivates  of  E^  and  H  , 
as  well  as  E  sni  H  themselves  have  finite  discontinuities.   If,  however,  the 
discontinuities  in  E„  PJid  H  have  this  form  bat  any  of  their  derivatives  become 

0       0 

infinite  at  a  discontinuity  tgj  of  t,  the  continued  integration  by  pairts  cannot  be 
carried  out.  Such  is  the  case,  for  exaniple,  when  the  form  of  B  or  H  near  t^ 
is  that  of  the  function  z  for  z  near  0  where  m  is  fractional.  Here  continued 
differentiation  of  E  or  H  will  ultimately  produce  infinite  discontinuities, 

0      0 

To  handle  these  more  complicated  discontinuities  a  new  form  of  the 
steady  state  field  will  be  introduced  which  will  likewise  give  rise  to  an 
asymptotic  expression  but  >diich  avoids  the  use  of  higher  derivatives.  The 
new  form  is  embodied  in  the  complex  integral  representation  (10.12).   It  will 
be  sliown  by  mathematical  examples  how  this  representation  can  be  used  to  obtain 
asymptotic  representations  of  steady  state  fields.  The  first  of  these  examples 
presents  a  discontinuity  which  can  be  handle<jl  by  the  method  of  the  preceding 
article;   the  others  cannot.  Among  these  examples  is  one  of  considerable 


-57- 


generality,  namely,  the  result  incorporated  in  the  series  (10. 3O).  The  type 
of  discontinuity  treated  "by  this  series  does  arise  in  an  application  to  diffrac- 
tion which  will  he  treated  in  a  later  paper. 

10.2  Let  us  assume  In  this  section  that  E  and  H  are  sectionally 

o      o     

analytic.  This  laeans  that  E  and  H  possess  continuous  derivatives  of  any 
0      0 

order  in  the  interior  of  any  domain,  D,  of  the  x,y,z,t  space  which  does  not 

contain  any  points  of  the  wave  front 

hyper  surfaces  "^^  -  ct  =  0  or  of  the 

STirfaces  of  refraction.  This  assumption 

ahout  E  and  H  is  Justified  if  fe  ^ /u   ,^ 

and  the  components  of  the  vector  field  F 

are  sectionally  analytic.  While  E  and 

Fig.  15  H  will  have  finite  jump  discontinuities 

0 

at  any  wave  front  or  surface  of  refraction,  we  do  not  assume  that  the  higher 

order  derivatives  of  E  and  H  attain  finite  limits  if  one  of  these  surfaces 

0      0 

is  approached. 

We  write  the  formulae  (9.15)  in  the  following  modified  form 


>X 


(10.1) 


u  = 


g  = 


lim 
T-^o 


lim 
T-*oo 


{e^(T)  e^^  -  i«f  E^(t)  e^'*'*  dt] 
fs^d)  e^^^  -  io/'  E^(t)  e^*^*  dt]  . 


We  next  define  two  vector  fields  A(x,y,z,t),  B(x,y,z,t)  "by  the  integrals 


A 


2ni  J 


00 


(10.2) 


B 


2ni  J 


00 


(r-t)2 


(T-t)' 


dT 


dT 


and  interpret  t  =  t.  •♦•  i  tp  as  a  complex  parameter.  Thus  A  and  B  are  complex 
functions  of  the  complex  variable  t  (and,  of  course,  of  x,y,z,  which  however 
are  fixed  for  present  purposes). 
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The  properties  of  these  complex  functions  are  test  studied  "by 

applying  a  formula  for  integration  by  parts  which  holds  for  Stieltjes  inte- 

grsJ.8  and  which  thereby  permits  us  to  treat  E  (7*)  even  though  its  derivate 

o 

may  have  infinite  discontinuities  on  the  range  0  to  oo  for  f  ,     According  to 

this  formula   A  and  B  may  be  written 

00   d  E  (  r) 
0 

r-  * 

(10.3)  ,  , 

B(t)  5^ 


A(t)^ 


o 


0 


r-  t 


where  the  integrals  are  Stieltjes    integrals.     For  t  =  t.   •♦■  it     not   on  the 

positive  real  axis,  A  and  B  are  analytic  functions  of  t.      If  E     and  H     are 

0     o 

sectionally  analytic  in  T*  the  analytic  functions  represented  by  the  two  inte- 
grals when  to  ^^  0  can  be  continued  past  the  positive  real  axis  along  paths 
which  cross  the  axis  at  values  of  t  where  E  (T)  sx^i^  H  (T)  are  analytic  in 
the  real  variable  T.  However  at  values  of  T  ,  where  E  and  H  are  discon- 

0       0 

tinuous,  the  functions  A(t)  and  B(t)  possess  singularities.  The  analytic 

continuation  through  values  of  "^where  £  and  H  are  sectionally  analytic 

will  createmiltiple-valued  functions  of  A  and  B  one  branch  of  which  in  the 

upper  half-plane  (t„  >  0)  is  represented  by  the  integral  and  another  branch 

of  which  is  represented  in  the  lov/er  half-plane  by  the  integral  when  X     <   0, 

Since  the  singularities  of  A(t)  and  B(t)  are  at  those  points  t^  -    —y^ 

on  the  positive  t.-axis  where  E  and  H  are  discontinuous,  the  stretches  of 
1  0      0  ' 

the  positive  t.-axis  between  singularities  may  be  taken  to  be  the  branch  cuts 
of  the  Riemann  surfaces  for  A  and  B.  Also,  since  the  negative  t.-axis  is  a 
region  of  analyticity  for  A(t)  and  B(t),  one  may  continue  from  t  >  0  to  t_  <  0 
across  the  negative  t,-axis  and  remain  on  the  same  branch  of  the  ftinction.  By 
following  such  a  path  one  passes  from  the  integral  representation  of  the  func- 
tions in  the  upper  half-plane  to  the  integral  representations  in  the  lower  half- 
plane  . 

11.  The  forimila  for  integration  by  parts  referred  to  is  proven  in  Wintner,  A: 
Spektraltheorie,  S.  Hirzel,  Leipzig,  1929.  The  assertions  which  follow  in 
this  paragraph  regarding  A  and  B  are  all  established  for  the  equivsilent 
forms  (10.3)  in  Wintner,  where  a3.so  a  full  discussion  of  the  Stieltjes 
integral  may  be  found.  Pages  83  to  88  and  pp.  9I-IOO  are  most  relevant 
background  material  for  this  paragraph.  Vintner's  integrals  are  defined 
over  the  range  -00  to  00.  Since  our  E  and  Hq  are  0  for  t  <  0  we  may 
restrict  the  integrals  here  to  the  range  0  to  00 .  Some  of  the  material  in 
Wintner  may  also  be  found  in  Widder,  D.V.:  The  Laplace  Transform,  Princeton 
University  Press,  I9UI,  p.  32U  ff. 
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One   can  verify  readily  that  the  vectors  A  and  B  satisfy  Maxwell's 
equations   {2.2),   namely, 

curlB-     -^A,     -     *L£a=     1    0, 
etc  c       t 

(10.1+) 

M 
curl  A  +     —  B.      =0 
c      t 

where  t  is  interpreted  as  a  complex  parameter  and  where  G  is  ty  definition  the 
vector  field 

since  for  the  pulse  solution  F(T)  =  g(x,y,z)/^  (T)  v;here  *7  (^)  is  given  "by  (9.1+). 
Indeed,   "by  differentiating  formula   (10.2)  under  the  integral  sign,    remembering 
that  Eq  and  H^  are  discontinuous  at  T=  ^    (cf.    the  derivation  of  formula  (8.U) 
we  obtain 


grad{^    4\To< 


'00       curl  E 


curl  A  =  -  3-r  >       ^i ^ ••■     ^5^  /  ^o     dr 

(10.6) 


and  similar  equations  for  curl  B,  j^-r    and    ^-r     . 

^  t       9  t 

The  equations  (10.1+)  thus  are  consequences  of  Maxwell's  equations 

for  3  ,  H  and  of  the  discontinuity  conditions  (3.9). 

10.3  It  is  possible  to  invert  the  integral  formulae  (10.2)  provided  that 
B  (t)  and  H  (t)  are  sectionally  continuous.   Let  f ,   as  "before,  denote  real  values. 
Then  according  to  Stieltjes'  inversion  theorem: 

e>0  Jo   (""^h   *  ^^^  -  ^^h  -  ^€  ))  dt^  =  E^(7-)  -  B^(0) 

(10.7) 

€-i^     K  C^^H  *  ^^^  -  ^^*1  -  ^^^^    **1  =  ^o^^^  -  ^0^°^  • 
12.   See,  Wintner,    loo.    cit. ,  p.   96  or  Widder,   loc.    cit. ,  p.    339. 
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If  E  or  H  are  discontinuous  at  the  point  7*  then  E  if)   and  H  (  7*)  are  to  be 

0       0  0  0 

replaced  "by  the  arithmetic  means 
(10. g) 


-  (E  (r+  0)  +  E  (r-  0) 

20  0 


I  (H^(r+  0)  -^H^cr-  0)    . 

Consider  now  the  first  of  equations  (10. 7).  ¥e  multiply  by 
iCOT 
iCOe  and  integrate  with  respect  to  T  from  T's  0  to  7*=  T.  By  carry- 

ing out  the  integration  under  the  limit  sign  and  integrating  by  parts  on  the 
left  and  integrating  directly  the  second  tenn  on  the  right  we  get 


€ 


i.<Jjf 


e      (A(r+  i£)-A(r-it))dr/ 


i«ilT 


=  ia>/   B^(r)  e"^'   dr  *  E^(0)  (1  -  e^-'^")  . 

0 

This  gireb  with  the  aid  of  (10. j)  and  (10. l): 


The  limit  on  the  left  side  is  equal  to  the  integral, y  A(t 

L 


N  i^t 

^  e    dt,  where  the 


path  of  integration,  L,  is  a  ctirve  enclosing  the  positive  real  axis  as  illustrated 
in  Fig.  l6  .  Since  E^  is  discontinuous  at  7* =  0  we  must  replace  E  (O)  in  (10.9) 

by  ^  E  (+  0).  The  latter  however  by  (3.11)  is  -  -^  "t(x,y,z).  We  get  there- 


2  o 
fore 


(10.10) 


2€ 


u  =  -  ^  i*  J  kit 

L 
V  =  J  B(t 


)   e     dt 


)  e     dt. 


i 


^  -toil  tst.x 


Pig.  16 

10. U  We  may  deform  the  path  of  integration  L  and  the  cuts  along  which 
the  branches  of  the  Riemannian  surface  are  joined  together.     The   integrals 
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(10,10)  will  not   change  if,   after  the   deformation,    the  path  L  still  surrounds 
all  the  csuts  of  the  surface."-^     This  leads  us  to  Fig.  f^.      The  cuts  go  from 


,.1 


v.y 


'    .i 


v_y 


B 


1 


^17 

rig.  17 


w 


each  singularity  t^  to  t(x    +  i  a> 
and  the  path  L  surrounds  all  these 
cuts.     Novf  consider  a  point  x,y,z 
in  whose  neighhorhood  g(x,y,z)  =  0. 
Then 


u  =  /  A(t)   e^^^^^     dt 


t, 


(10.11) 


■I« 


t)  e  dt  . 


Since  at  such  an  x,y,z  Eq,  H  and  all  time  derivatives  are  zero 
at  t  =  0,  in  this  case  we  conclude  that  there  is  no  singularity  at  t  =  0, 

We  now  contract  the  different  sections  of  the  path  L  to  the 
straight  lines  which  represent  the  cuts.  We  replace  the  variahle  of  integra- 
tion t  hy  —  (1/    •*■  T  and  obtain  u  and  v  in  the  form 


(10.12) 


where 


ot=l 


Ve<  e 


ikf^ 


M 


T«  -.  [  ^(\^fi,  .r)  .^r^r. 


M 

where  the  path  M  of  integration  now  surrounds  the  positive  imaginary  suxis  of 
the  complex?'  -  plane. 


13,      Strictly  the  path  L  should  first  he  deformed  into  finite  U  shaped  paths 
joined  by  Jorizontal  sections   such  as  AB  in  the  figure  and  then  let  AB 
approach  infinity.     The  function  A(t)e^**^   approaches   zero  however  as 
t  — »  a  +  ioo  for  any  a.      This   follows  from  the  fact  that  A(t)  which  is 
given  by   (10.2)   approaches   0   as   the  imaginary  part  of  t  >    i    00     and  so 

does  8^**'*    •     EeiiceJ      A(t)  e^*"*  dt  approach  aero  as  AB  moves  upv/ard. 


AB 


w 
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The  result   (10.12)  is 
the  "basic  result  of  this  article 
^ z  T  "^  ^'fi.  ^^'^  famishes  the  second  method. 

of  ohtaining  the  asymptotic  ex- 

M 

pans  ion  of  the  'steady  state 


%  field.  The  actual  form  of  the 

^^S*  ^^  result  does  not  make  evident 

two  facts  ahout  the  expansion  more  evident  in  the  formulation  of  article  9,  The 
values  of  the  Ug^   and  Vo(   ,  since  they  are  complex  integrals  over  closed  paths, 
will  actually  depend  upon  the  singularities  of  the  integrands  A  and  B.  However 
it  was  pointed  out  earlier  that  the  singularities  of  these  fujictions  are  those 
of  the  pulse  solution  I  and  H  .  Hence  the  pulse  solution  is  involved  as  in  the 
earlier  formulation.   The  other  essential  feature  of  the  previous  asymptotic  ex- 
pansion is  that  the  hasic  variable  is  the  wave  length  X.     It  will  soon  appear  that 
each  U^  sad.  ^  in  formula  (10.12)  is  an  infinite  series  in  ^  or  in^so  that  again 
the  final  expressions  for  u  and  v  are  asymptotic  expansions  in  A. 

10.5  ^iie  individual  terms  of  the  series  (10,12),  namely, 

(10.13)  ^^  =Uo(  e^^^ 

are  solutions  of  Maxwell's  equations.   Indeed 

curl  Vol  =  J  (curl  B  +  -  grad^  J<     ■^)  e^^'^'dr 


M  -  ■  9*' 


=  / 


curl  Be  d^  -  ik  gradlZg,X  le(  . 

M 

Hence  "by  (10.1<-)  and  partial  integration,   G  "being  0  at   (x,y,z), 
curl  V^     +  ik  grad|4i  ^  V«      =     — CrU„     -  ik£  Uoc     , 


c 


(10JJ+) 


and  hence 


(10.15) 


curl  Vec    +  ik  grad^^X   \Joc    =     ik/<Vac 
curl  Tfej     ■•■  ikfiuoc      -     ^^CT-a^    =  0, 
curlOot      -     ikyU  -^  =0. 
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10,6  In  the  next  few  pages  we  shall  illustrate/the  general  result  (10.12) 
can  "be  applied.  The  illustrations  are,  in  so  far  as  this  paper  is  concerned,  with- 
out physical  significance.  However,  in  a  paper  on  diffraction  now  "being  prepared, 
it  will  appear  that  the  second  of  the  illustrations  pertains  to  physical  prohlems. 

The  first  illustration  deals  with  an  E  and  H  whose  discontinuities  and  the  dis- 

0      o 

continuities  of  their  derivatives  are  finite  jumps  and  could  therefore  he 

treated  by  the  method  of  article  9»   Ti^e  second  illustration  nrast  necessarily 

he  treated  hy  the  method  of  this  article. 

Consider,  for  example,  the  case  that  B  and  H  and  all  their  deriva- 
'  ^       '  0       0^ 

tives  approach  finite  limits  if  "T  approaches  0,  that  is,  if  t  =  —  ^^(x,y,z)  ■♦•  0, 

Let  us  also  assume  that  E  =  H  =  0  for  t  =  —  ^f^  -  0  so  that  t  =  —  If^   represents 
the  directly  transmitted  wave  front  of  the  pulse  solution,   We  may  expand  E  (t)  in 
the  neighborhood  of  t  =  —  tt^  by  means  of  Taylor's  theorem  for  t  values  greater 

than  ^—=-    ,  Then 
c 


:«)=!  Si^-lr^.r  .  w^ 


(10.16)  V''   piO  tTT  --c'-i'   ■     '■    'c 

wherein  E  means     ^ 


^att-^l^,tO. 


For  any  E^Ct)  the  corresponding  A(*t  )   is  given  by  (10.2).     Now  the 
behavior  of  the  function  A(t),   in  the  neighborhood  of  T"-  0  or  t  =  "Z^i     ^^ 

the  same  as  that   of  the  integral 

_     -..  ._       .  .  4^ 

III/:   i/«o 

where  w   is  arbitrarily  small.   Since  we  are  only  interested  in  the  singular  part 
of  A(t)  v;e  are  permitted  to  omit  regular  contributions  of  the  integral.   To  deter- 
mine the  singular  part  of  expression  (10, I7),  replace  t  -  f\)vfit   say,  and  reverse 
summation  and  integration.  The  integral  then  has  the  form 


tne  same  as  tnat  01  ine  ini^egraj. 

UO.^7)      A,U,=  fe^J,,,,  f^        (T-  » 


I 


— 2 —  r 
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Since»/ls   integral,   'binomial  expansion  of  the  numerator  shov/s   that  the  only 

singular  term  will  "be  — .      We  can  integrate  this   singular  term,    for 

the  parentheses  are  independent  of  JQ  ,  Taking,    then,    the  factor  (t  -     ^r    1)     out 
of  the   integral  we  may  express   the  singularity  of  the  integral  at   t = Yy c  merely 
by  log(t  -yi/c). 

This  result  leads  to  the  following  expression  for  the  singularity 
of  A.  (t)   or  of  A(t)    : 

*^ 

We  may  introduce  this  expression  into  the  "basic  result  (10,12)  with- 
out changing  the  result  for  Uof  hecause  the  integral  of  this  regular  part  around 
the  closed  path  M  of  Figure  /^will  be  zero.  It   follows  since  f=  t  -  r'lc 

that  v—      I         r/ 

^««    =  2fu  u=o  ^^  vf-*    log  r  e        dr . 

.  Since  log  T  represents  the  singularity  of  the  term  at  7*=  0  or  at 

t  =    ^- —    and  since  the  values  of  log 7'  on  the  two     sides  of  the  branch  line 
c 

differ  by  the  constant  -2nl  we  conclude 

■loo       E 


CO  -.  -  -V 


(10.19) 


-TX.       /•loo        E 


Uec    =    -ili>2_    /  ^T^'e^^dr 


-Z:    Jrlf   ,^'' 
*'=°      (ia>r    ° 


The  same  result  follows   dlrecly  from  (9.18).      It   should  be  noted 
that  E     =  /  B    I  since  the  left  hand  derivative  is     0  In  this  example. 

/I  ftf^^ 
F(^)  e  d^  where  FC^')  is  regular  is  zero  since  the 

M 

path  is  from  1  oo  to  0  and  then  from  0  to  1  oo  . 

15.  Integrate  -y^'e^***  df   by  parts  letting  u  =  7'^    and  dV  =  e^**'''dT  . 

Repeat  this  V  times;   the  explicit  'hiv"  term  drops  out  each  time  leaving  finally 
the  integral  of  e-^*"  with  a  constant  multiplier. 
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We  now  illustrate  the  application  of  the  "basic  result  (10.12) 
with  an  example  in  which  the  method  of  the  preceding  article  does  not  apply. 
We  consider  the  case  where 

2 


c 

Again  A(t)  is  given  hy  the  integral  in  (10,2),  the  T'of  that  fonnula  being  equal 

to  t ^  in  the  representation  of  E  here.  The  hehavior  of  A(t)  in  the  neigh- 

1  i 
horhood  of  t  =  —  w^is  given  "by  the  integral, 

0  T*-  (t  -  -Mr) 

where  t  here  assumes  veilues  greater  than  "^  .  By  repetition  of  the  argument 

used  in  the  preceding  examples  to  obtain  formula  (10.18)  we  may  state  that  the 
singular  part  of  A-,(t)  is  given  by 

CO  if 

(10.20)  A,(t)  =  1^  ^  ^^  cjt  -  \'^r  J    ^J^[\^^      ^r. 


Jvr 


0 

c 


Since  f dT   has  no  singularity  at  7^=  0,  it  follows  that  the  singularity  of 


r 


the  integral  in  A,  (t)  is  the  same  as  the  singularity  of  the  integral 

( -L-- -  -l-jdT. 

0     r-(t  -  ij^)         Vr 

By  adding  fractions  and  by  repeating  the  reasoning  above  in  con- 
nection with  formula  (10.18)  we  have  that  the  singularity  in  this  integral  is 
the  same  as  the  singularity  in 

(t .  ^)  Y l£ 


-6&- 

By  adding  and  subtracting  the  integral  for  the  range  O  to  oo  we 


have  that  this  last  integral  equals 


(t  -  -*i_)  I        ^ rTT"^  *  Regular  terms. 


c 

The  integral  can  he  evaluated  hy  a  standard  integral  form  and  we  ohtain 

=  -  TTi  Vt  -  Je^   +  Regular  terms. 
c 

By  substituting  this  result  in  formula  (10.20)  we  obtain  that 


A(t)  =  -  i  %-    ZZ   C  (t  -  2^f*'^/^    +  Regular  terms. 


2  <i*  *^0  ^ 
We  now  introduce  this  value  in  (10,12)  and  obtain  for  U^ 

.n+1/2    iCuT 


-  7"  c     f    i  i- 


n=0   ^  M    "^  ^'^ 

Now  the  path  M  goes  from  i  »8  to  0  on  the  right  side  of  the  branch  line  and 
from  0  to  iOO  on  the  left  side.  Because  the  quantity  ^       '      =  r  7"  '  "^ 
appears  in  the  integrand,  the  integrand  is  multiple-valued  by  virtue  of  the 
7*  '2  .  On  the  tv^o  sides  of  the  branch  line  the  two  values  ot  T   '      differ 
by  a  negative  sign.  Hence  we  may  combine  the  two  integrals  thus: 

^       n=0   ^•'o      ^ 
By  a  change  of  variable  we  obtain 

^    (-ic«;)"^"'"«'o 


An  integration  by  parts  gives 


2       F    — 

e-^   dr 
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In  view  of  the  definition  of  the  Gamma  function  we  get 

(10.21)  ^=2.    "^^TTp        JL  (^*3/2)  . 

This  result  shows,  as  did  our  preceding  example  that  Uot  is  an  infinite  series 

in  the  hasic  vajriahle  —  or  in  ^  eind  it  shows  farther  that  Ugn  approaches  0  with  X, 

10,7  Out  examples  show  how  the  singularities  of  A(t)  and  B(t)  determine 
the  U^  and  Vpc  and  hence  the  u  and  v.  Since  the  analysis  is  lengthy  it  is 
desirahle  to  obtain  a  general  result  on  the  value  of  U<x  and  "'Igi    for  a  general 
class  of  functions  A(t)  and  B(t)  and  thereby  avoid  frequent  recourse  to  such 
analysis.  Me   therefore  now  shall  obtain  a  general  result  of  this  nature.  It  is 
this  general  result  which  will  be  applied  in  a  later  paper. 

Let  us  assume  that  Sp(t)  and  H  (t)  are  sectionelly  analytic.  Consider 
for  example  E-Ct)  in  the  neighborhood  of  a  point  t  of  discontinuity.  Let 


"^  C 


E^  =  CHit     -  t)  2_  -^^K-  ^^"^       .  for  t  <  t^ 
o      0      ^-  ,   n.'   0        '         0 
n=l 

(10,22)  e>o     ^ 

%-0^''-'o^    §1  ST^*-  ^o^""    .  fort  >t^  , 

whereOt  and  4  are  functions  of  t  -  t.  We  shall  require  that  the  functions<V('7') 

sxid.pif')   be  analytic  functions  of  /'which  are  not  regular  at  7'=  0,  for  example 

k 
algebraic  functions  of  the  type/  ,  where  k  is  not  necessarily  an  integer.  We 

also  include  the  case  that  0{  or^tf  equals  log  7^  . 

We  knov;  that  A(t)  as  defined  by  (10.2)  has  the  same  singularities  at 


t  =  t  as  the  sum  of  the  integrals 


■^  V 


,    /"'o     Edr  ,     /■to'' 

whereQ  is  arbitrarily  small. 


-68- 


We  consider  the  terms  of  the  series  (10.22)  individually,  i.e.  ve 
consider  the  functionOC(t  -  t)(t  -t)'^  for  t  <  t  and^(t  -  t  )(t  -  t  )'^  for 

0         0  f    ^  0         0 

t  >  t  .   The  contrihution  of  these  terms  to  the  integrals  in  (10,23)  are 
o 


^"'"  =  -  h  ^.f-^Hf. 


0  r-(t^-t) 


^     r^  dr 


V^''   dtef^J  r-(t  - 1  )  ^    <^^- 

0  0 


By  the   same  argument  as   that  used  to  derive   (10,18)  from  (10, 17)  we  conclude 


that  ^Y  and  A     differ  only  "by  regular  terms   from  the  expressions 


0     '       0 

(10,?1|)  f^     JO   f** 


0  •      o 
We  next  determine  the  contrihution  of  these  expressions  to  the 


function 


(10.25)  U=  /    A(t^  +70  e^^^dr  . 

Here  U  denotes  any  one  of  the  ^ki  in  (10,12).     Letting 


(10.26) 


'0 
0    '  0 


we  get  from  (10,26),    (10. 2U)  and  (10. 25) 

U^  =  -ia;/    f(t)t^  e"^***^     dt 


(10.27) 


/ 


u;=-iCuJ    g(t)t^e^«'^     dt 
"^  M 
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'  Fig.  19 
By  introducing  the  differences 

[f]=  fi  -  fj 


where  M  and  M'  are  the  paths  of  integra- 
tion indicated  in  Fig./3«   These  two 
paths  surround  the  positive  and  negative 
imaginary  axis  respectively.  The  two 
"brajiches  of  each  path  lie  in  different 
sheets  of  the  Riemannian  surfaces  of  f  and 

g* 


fg>  gl  - 


Sr 


in  corresponding  points  of  the  two  paths  (Fig./9)  we  may  transform  (10. 27)  into 

'-i  00 


u;;  =  -  icu 


(lo.zg) 


r  f  0 


^n  -i«t 
t  e 


dt 


'0 
ri  00 


example 


We  determine  the  differences  [f]   and  [gj   as  follows.   Consider  for 


X 


The  integral  determines  the  value  f^(t) 
in  the  first  sheet  of  the  Riemannian  sur- 
face. The  associated  value  f2('t)  ^^  ^^^ 
second  sheet  is  obtained  Toy   the  integral 


f( 


^^  ~     2ni   <       7"-  t 


Fig.    20 


v;here  X  is  a  curve  connecting  the  points 


0  andO   such  that  the  point  t  lies  hetween  the  curve  ^  and  the  real  axis.     By 
Cauchy's   integral  theorem  we  conclude 


ff]   =  fj^Ct)    -  f^Ct)  =0C(t) 
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and  similarly 

g  =  g^Ct)  -  g2(t)  =^(t)  . 

We  introduce  these  restilts  in  (10.28)  and  obtain  after  changing 
the  variahles  of  integration: 

'        .n       -t     ^^ 
t       e         dt 


"^-it^jt*^'^' 


(10.29) 


Since  u.  and  u  are  the  contriliutions  to  U  of  A?  and  A^  ,  the  con- 

trihution  of  A  itself  to  U  or  the  complete  asymptotic  series  of  the  foinction  U 
thus  has  the  form 

if  E^,  in  the  neighborhood  of  t  =  t^  =  c^'  ^^  given  ty  the  expressions  (10.22). 

Of  course  the  U  so  obtained  is  typical  of  each  of  the  tW  to  be  evaluated  for  the 
representation  (10.12)  of  the  steady  state  field. 

Consider  as  aa  example  the  case  that  0C=  t  '  p~  t       where  k  and^ 
are  greater  than  -1.  We  obtain  readily  by  use  of  the  Gamma  function  that 

D„J*(n+l+Cr)     C  r(n+l+k) 


(10.31) 


_  5"  1  f^Ti^+l-'tr)  C^r(n+l+k)\ 

""^  "^0   ^l  (-ia»)^*<^   ■  U^IiP^   )      • 


Consider  next  the  example 
C^  =  (-1)^  D^  ande(=y^  =  1  log  t.  This  meajis  that 


V')  =  iloglt-t„\r-i,Ct-t^)° 


on  both  sides  of  t  =  t  .  In  this  case  we  obtain  from  (10. 3O)  that 

(10.32)  TJ  =  i  ^-.  — ^  . 

n=0  {-xCof 

i.e.,  except  for  the  factor  i,  the  same  result  as  in  case 


-o=-^;iT(t-t^)^  ,      t  >  t^ 


=  0  .   t  <  t. 
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It  Is  apparent  from  these  examples  that  the  result  embodied  in  formula  (10. 30) 
can  "be   used  to  obtain  explicit  asjTnptotic  expansions  in  a  variety  of  cases. 
Further,  since  U  in  that  formula  is  typical  of  the  Uq^  In  (10.12)  it  follows 
that  we  can  obtain  asymptotic  expressions  for  the  steady  state  fields  u  and  v. 
To  obtain  terms  in  the  expansions  of  u,  say,  of  a  given  order  in  •; —  it  would 
be  necessary  to  collect  these  terms  from  the  various  expressions  for  the  Uoc 
Formula  (10,30)  will  prove  useful  in  applications  thus  far  made  of  the  above 
theory.   One  of  these  applications,  namely,  to  the  asymptotic  evaluation  of 
diffraction  integrals  should  be  completed  shortly. 


